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ABSTRACT 


Methods  are  explored  for  estimating  Range  and  Range-Rate 
from  the  geometry  and  kinematics  of  the  air-air  homing-missile 
combat  situation  rather  than  from  direct  measurements.  The 
methods  reported  herein  use  the  various  signals  available  from 
the  missile  autopilot,  plus  the  inertial  rotation-rate  of  the 
Line  of  Sight  provided  by  the  seeker. 

It  is  shov.T.  that  acceleration  of  the  target  in  the  radial 
direction, parallel  to  the  Line  of  Sight,  precludes  successful 
estimation  by  the  classical  techniques  such  as  adaptive  parame¬ 
ter  Identification,  Kalman  filter  estimation,  or  the  various 
minimum-variance  estimation  methods,  A  nonlinear  estimator  is 
described  which  estimates  the  target  acceleration  and  velocity 
components  and  can  yiela  accurate  range  and  range-rate  esti¬ 
mates  if  correctly  initialized  at  launch.  This  estimator  is 
able  to  take  advantage  of  the  many  inequalities  which  constrain 
the  maneuvering  of  an  air-air  target.  Errors  and  performance 
of  this  estimator  are  demonstrated  by  computer  simulations. 
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1.  INTRODUCTION 

Proper  use  of  range  and  range-rate  data,  or  the 
equivalent  "time-to-go, "  can  Improve  the  performance  of 
an  air-air  missile.  At  the  same  time,  instruments  whi-".; 
directly  measure  range  and/or  range-rate,  such  as  radars, 
are  expensive.  The  purpose  of  this  study  was  to  determine 
whether  adaptive  parameter  identification  techniques  could 
enable  estimation  of  range  and  range-rate  from  less  expen¬ 
sive  sensors  which  are  used  for  example  in  infrared 
seekers . 

1.1  Background 

The  trajectory  of  ar.  air-air  homing  missile  may  be 
regarded  as  the  response  of  a  closed  loop  to  an  external 
command.  In  this  case,  the  external  command  is  the  tra¬ 
jectory  of  the  target,  which  must  be  assumed  to  be  maneu¬ 
vering  either  to  conduct  its  attack  or  to  attempt  to  evade 
the  missile.  The  closed  loop  consists  of  the  geometry 
which  yields  the  Inertial  Line  of  Sight  (LOS)  rotation-rate, 
the  navigation  or  guidance  computer  which  forms  commands  to 
the  missile  as  a  function  of  the  LOS  rate,  and  the  missile's 
dynamic  response  to  that  command,  which  closes  the  loop. 

The  inertial  LOS  rate  Is,  as  viewed  from  the  missile,  a 
dynamic  function  of  the  target  and  missile  accelerations 
perpendicular  to  the  LOS,  and  the  range  and  range-rate. 
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If  signals  of  the  LOS  rate  and  the  missile  acceleration 
perpendicular  to  the  LOS  are  available  measurements,  it 
might  be  possible  to  deduce  the  dynamics  of  the  geometry, 
which  may  be  regarded  as  a  "plant"  within  the  closed  loop  of 
the  entire  system. 

As  stated,  the  purpose  of  this  research  study  was  to 
explore  means  of  estimating  the  missile-target  range  and 
range-rate  by  applying  adaptive  parameter  identification 
techniques  to  this  "plant"  or  geometry. 

This  approach  to  estimation  of  range  and  range-rate 
uses  the  missile  dynamics  as  a  probe  v.’lth  which  to  examine 
the  geometry  or  kinematics  of  the  system.  It  is  therefore 
called  "Dynamic  Ranging,"  or  "Kinematic  Ranging." 

Figure  1.1,  right,  shows  the  closed  loop  from  target 
acceleration  through  the  kinematics  and  dynamics  to  the 
missile  response,  In  the  upper  portion.  The  lower  part  of 
the  diagram  shows,  in  a  symbolic  form,  the  parameter  esti¬ 
mation  computer  whose  ultimate  outputs  are  the  target  range 
and  range-rate.  Consider  the  "box"  containing  the  geometry 
and  kinematics;  its  output,  the  LOS  rate,  is  a  measurable 
signal,  while  the  two  Inputs  are  the  target  and  the  missile 
accelerations  perpendicular  to  the  LOS,  But  while  the 
missile  acceleration  is  an  available  signal,  the  target 
acceleration  is  not;  estimation  of  range  and  range-rate 
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therefore  must  depend  significantly  on  the  characteris¬ 
tics  of  the  target  acceleration  perpendicular  to  the  LOS. 

It  will  he  seen  that  the  target  acceleration  parallel  to 
the  LOS  also  has  a  strong  effect  on  the  estimation  problem. 
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1.2  Organization  of  the  Report 

The  principal  results  of  this  study  and  the 
conclusions  are  gathered  in  Section  2.  The  problem  is 
described  and  mathematically  posed  in  Section  3.  A 
nonlinear  range  and  range-rate  estimator  is  described 
in  Section  4,  together  with  an  error  and  performance 
analysis.  Several  linear  estimators  are  described  in 
Section  5>  together  with  comments  on  the  conditions  for 
which  the  problem  is  mathematically  observable.  A 
technological  forecast  is  presented  as  Section  6, 
References  are  gathered  in  Section  ?•  Appendix  A 
contains  some  details  of  the  derivation  of  the  nonlinea 
estimator,  and  Appendix  B  contains  tabulations  of  the 
time-varying  coefficients  of  its  error  propagation  equa 
tion  for  a  variety  of  typical  trajectories. 
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SUMMARY  AND  CONCLUSIONS 

A  summary  of  the  findings  of  this  study  is  presented. 
The  principal  conclusions  are  stated  and  discussed. 

2 . 1  Summary 

Methods  are  explored  for  estimating  Range  and  Range- 
Rate  from  the  geometry  and  kinematics  of  the  air-air 
homing-missile  combat  situation  rather  than  from  direct 
measurements.  The  methods  reported  herein  use  the  various 
signals  available  from  the  missile  autopilot,  plus  the 
inertial  rotation  rate  of  the  Line  of  Sight  provided  by  the 
seeker. 

It  is  shown  that  acceleration  of  the  target  in  the 
radial  dlrectlon^parallel  to  the  LOS,  precludes  successful 
estimation  by  the  classical  techniques  such  as  adaptive 
parameter  identification,  Kalman  filter  estimation,  or  the 
various  minimum-variance  estimation  methods ,  unless  the 
sensor  noise  levels  are  very  low  and  the  geometry  of  the 
combat  is  favorable.  A  nonlinear  estimator  is  described 
which  estimates  the  target's  aspect  angle,  its  velocity  and 
acceleration  components  plus  turn-rate  and  roll-rate, and 
can  yield  accurate  range  and  range-rate  estimates  if 
correctly  initialized  at  launch.  This  estimator  is  able  to 
take  advantage  of  the  many  inequalities  which  constrain  the 
maneuvering  of  an  air-air  target.  Errors  and  performance 
cf  this  estimator  are  demonstrated  by  computer  simulations. 
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2.  a  Conclusions 

1.  Target  trajectories  which  have  significant 
acceleration  parallel  to  the  line  of  sight  occur  in  most 
combat  geometries  except  head-on  and  tail-chase  configu¬ 
rations.  If  the  target  has  significant  acceleration 
parallel  to  the  line  of  sight,  then  range  and  range-rate 
are  mathematically  unobservable  from  the  kinematics  by 
linear  methods. 

In  practice  this  means  that  range  and  range-rate 
can  be  estimated  only  in  favorable  geometries  and  while 
active  maneuvering  of  the  missile  occurs  and  that  the 
sensor  noise  levels  must  be  so  low  that  very  short  aver¬ 
aging  times  in  the  estimating  filter  will  yield  acceptable 
accuracy.  Filter  averaging  times  of  the  order  of  one- 
quarter  second  to  one  second  yield  range-rate  errors  as 
large  as  30  to  50?  when  the  target  has  3  to  6g  acceleration 
parallel  to  the  LOS. 

2.  Two  approaches  are  available  to  resolve  these 
difficulties.  A  method  which  uses  linear  principles 
combines  the  kinematic  techniques  with  stadlmetric  methods. 
Stadimetric  methods  require  that  the  seeker  be  able  to 
measure  some  function  of  the  target's  area  or  angular  size 

and  rely  on  the  assumption  that  the  target  is  of  constant 
■  tze.  Using  stadimetric  data  relieves  but  does  not  eliminate 

the  difficulties  noted  in  conclusion  (1),  above.  A  second 
method  consists  of  a  nonlinear  estimator,  discussed  below. 
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3.  A  nonlinear  estimator  which  estimates  the  target 
velocity  vector  components  and  acceleration  components 
relative  to  the  line  of  sight  and  also  estimates  the 
target  turn-rate,  and  its  rate  of  change,  can  be  devised. 
As  a  stand-alone  unit,  this  estimator  requires  Initiali¬ 
zation  of  range  end  range-rate.  In  combination  with  other 
equipment  it  can  be  used  as  a  filter  to  exclude  impossible 
target  accelerations  due  to  seeker  noise  or  noise  origina¬ 
ting  in  some  other  method  of  ranging.  Its  principle  of 
operation  is  baaed  on  the  constraints  on  the  target 
maneuverability  dynamics,  and,  in  particular,  on  the 
fact  that  the  target  airspeed  cannot  be  significantly 
changed  during  an  alr-alr  missile  engagement. 
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3.  THE  MATHEMATICAL  MODEL  OF  KINEMATIC  RANGING 

The  study  is  restricted  to  the  Maneuver  plane. 

The  key  assumption  guiding  and  restricting  this  study 
is  that  the  seeker  in  the  missile  provides  the  line-of- 
slght  (LOS)  rate  only;  it  does  not  yield  any  indication 
of  target-size,  or  shape  or  changes  thereof,  as  television 
or  mosaic-type  infrared  seekers  can,  nor  does  it  provide 
any  range  or  range-rate  information  directly,  as  a  radar 
does,  nor  even  yield  the  variations  of  intensity  of  signal 
from  the  target.  The  simplest  type  of  seeker,  simply 
pointing  at  the  target  and  carrying  a  rate  gyro  to  sense 
the  LOS  rate,  is  assumed.  On  the  other  hand,  the  missile 
is  assumed  to  have  a  pair  of  accelerometers  in  the  maneuver 
plane,  so  that  its  accelerations  may  be  resolved  about  the 
line  of  sight.  Similarly,  it  is  assumed  that  missile  air¬ 
speed,  or  its  estimate,  can  be  provided  if  needed  and  that  no 
data  are  transferred  from  the  launch-airplane  to  the  missile 
after  the  missile  has  been  launched.  The  two  areas  which 
describe  the  problem  are; 

(a)  The  geometry  and  the  equations  which 
define  the  kinematic  relations  from  which  we 
hope  to  estimate  Range  and  Range-Rate,  and 

(b)  The  constraints  on  the  target. 
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These  areas  are  discussed  in  detail  below.  The 
problem  v/hich  we  attack  is  estimation  of  Range  and  Range-Rate 
from  the  kinematics  and  dynamics  of  the  combat. 

3.1  Kinematics 

Under  these  circumstances,  a  polar  coordinate  geometry 
with  the  moving  origin  located  at  the  missile  is  the  logical 
coordinate  system.  Figure  3-1  shows  the  geometry. 


Figure  3*1  -  Geometry 


We  first  derive  expressions  for  the  velocity  components 
parallel  and  perpendicular  to  the  LOS.  Similar  expressions 
for  the  accelerations  follow. 


The  instantaneous  linear  components  of  velocity  parallel 
to  the  Line-of-Slght  (LOS)  yield  the  polar  coordinate  com¬ 
ponent  , 

f  *  cos  (Y.p-o)-Vj^^  cos  (Yy,-'’) 


(3.1 
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while  the  components  perpendicular  to  the  LOS  are 
r  o  »  sin  (Y,p  " 

(3-2) 

The  angle  may  be  called  "The  Target  Aspect  Angle" 
as  it  represents  the  attitude  of  the  target  velocity 
vector  relative  to  the  LOS,  and  is  closely  related  to 
the  target  body  attitude  with  respect  to  the  LOS. 

As  missile  acceleration  components  are  assumed  to  be 
measured,  we  differentiate  (3-1)  and  (3.2)  so  that  the 
acceleration  terms  will  become  expllcity  evident, 
yielding 

r  =  CV^  cos  (Y,p-a)  -V^  sin  ) ]-[ cos(Yj^^.o) 

sin  J+c[V^  sin  sin  (Yp,,-CT )  ] 

(3.3) 

and 

ro+ro  =  [V^  sin  ( Yrp-o )+V„y,j,  cos  (Y.p-o)3- 

o[V^,  cos 

(3.^) 

The  components  of  acceleration  of  target  and  missile 
parallel  to  the  instantaneous  llne-of-sight  are  the  first 
two  terms  of  the  right  of  (3.3),  respectively. 
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Therefore  define  cos  (y.j.-o )-V^Y,p  sin  (y^-o) 

^  (3.5) 

r 

(3.6) 

The  third  term  on  the  right  of  (3*3)  represents  the 
apparent  acceleration  due  to  the  LOS  rotation-rate. 

Now  multiply  (3.2)  by  o  and  substitute  (rb  )  for  the 
third  term  on  the  right  of  (3.3);  with  this  substitution 
and  the  definitions  (3.5)  and  (3.6),  (3.3)  yields 


r  -  ro'  -  a^  +  a,^  =  0  . 

r  r 

Equation  (3.^)  is  similarly  simplified. 

As  before,  define  the  components  of  target  and  missile 
accelerations  across  the  instantaneous  line-of-slght  by 


^T^“Vt  sin  (Y.j,-o)+V^Y,p  cos  (Y,p-o) 

aj^  -Vj^  sin  (y^-o)+V^^;j^  cos  (y^^-o) 
0 


(3.8) 

(3.9) 


and  substitute  these  for  the  first  and  second  terms  on 
the  right  of  (3.^).  Again,  the  thli’d  term  on  the  right 
of  (3.^)  is  the  apparent  acceleration  due  to  the  LOS 
rotation-rate.  Now  multiply  (3*1)  by  o  and  substitute 
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(f  o)  for  the  third  term  on  the  right  of  (3.^),  yielding 
the  acceleration  equation 

r  e  f  2  ^  i  =■  0. 

a  o 

If  we  use  the  notation 


then  the  velocity  equations  (3.1)  and  (3.2)  may  be 
expressed  as 

r  -  cos  +  V|^  cos  =  0 
r  d  -  \.p  sin  +  Vj^  sin  =  0 

and,  with  the  more  compact  notation 

r  r 

’ 

0  o 

for  the  components  of  target  and  missile  velocity 
parallel  and  perpendicular  to  the  LOS,  we  have  the 
velocity  equations  in  the  form 


r  o  -  V„  +  Vj^  =  0  . 
0  'o 


(3.11) 


(3.12) 

(3.13) 


(3.1^) 


(3.15) 

(3.16) 
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Equations  (3.1)  and  (3*2)  or  their  equivalents  (3.12)  and 
(3.13)  or  (3.15)  and  (3.16),  together  with  (3.7)  and  (3.10) 
form  the  key  elements  of  the  mathematical  model  of  the 
kinematics.  The  assumed  noise  levels,  bias,  initialization 
errors,  etc.,  which  determine  the  errors  of  the  proposed 
solution,  are  gathered  in  Section 
•  * 

The  term  ro  in  (3-7)  deserves  a  brief  discussion. 

This  term,  a  component  of  the  range-acceleration,  is  due 

purely  to  the  geometry.  Its  physical  significance  can 

easily  be  perceived.  Assume  you,  (M),  are  standing 

motionless  at  the  coordinates  X  =  0 ,  Y  =  Yq  where  lY^j  is 

small,  so  that  =  Vj^,  *  =  aj,.  =  0.  A  vehicle,  (T),  is 

approaching  at  constant  velocity  on  the  X-axis;  X  may  be 

positive  or  negative  but  r  is  negative  as  the  vehicle  is 

approaching.  The  target  is  not  accelerating, so  that 

a^  “  start  of  the  problem,  r  ~  -  1X1, 

r  o 

while  at  the  end  r  Z  +  |x|;  Kange-Rate ,  f ,  has  changed 

from  -lx|  to  +  1x1;  the  polar-coordinate  acceleration  sigr.al 

.  2 

which  expresses  this  change  is  r o  .  If  the  offset  is 
very  sm;all,  the  change  from  -1x1  to  +1.\1  can  occur  very 
suddenly.  Implying  large  values  of  r.  In  an  air-air 
combat  situation.  It  is  therefore  possible  for  very  large 
range-accelerations  to  occur,  due  to  the  geometry  alo.ne, 
even  though  neither  missile  nor  target  maneuvers. 
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3.2  Constraints  on  the  Target 

The  constraints  under  which  the  target  airplane  operates 
form  an  essential  part  of  the  system.  An  airplane  can  increase 
or  decrease  its  airspeed  by  changing  its  thrust  or  by  use  of 
the  earth  gravitation  vector  component;  in  either  case  or  in  com¬ 
bination  the  longitudinal  acceleration  increment  is  not  more 
than  Ig.  An  airplane  can  also  change  Its  airspeed  by  increasing 
Its  drag;  this  can  be  accomplished  deliberately  by  use  of  speed 
brakes  or  involuntarily  as  a  result  of  the  Induced  drag  due  to 
the  lift  resulting  from  the  airplane's  evasive  maneuvers.  But 
the  target  of  an  air-air  missile  attack  must  try  to  maintain  air¬ 
speed  to  conserve  maneuverability.  We  therefore  consider  that 


is  a  plausible  constraint  for  deliberate  airspeed  changes,  and 
we  neglect  the  Involuntary  changes  in  this  study. 

The  target  turn-rate,  is  governed  by  several  inequality 

constraints.  At  high  airspeed,  jV^Y^I  is  limited  by  the  ability 
of  the  pilot  (or  structure)  to  withstand  g-loads,  equal  to 
VmTrp,  v;hlle  at  low  airspeed  the  turn-rate  is  limited  by  aero¬ 
dynamic  stall  related  to  bank  angle  In  addition,  the  target 

pilot  and  aircraft  roll-dynamics  establish  a  minimum  time  in 
which  a  maneuver  can  be  performed  or  changed  by  the  target  so  that 
I  YrpI  ,  or  1  ^,pl  ,  is  limited. 

The  high-speed  condition  provides  the  acceleration 
constraint,  due  to  pilot  or  structural  strength  limits, 
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Further,  turning  is  achieved  by  bank-angle.  Taking 
the  pilot's  characteristics  Into  account,  the  bank  angle,  4>rp, 
cannot  be  changed  from  hard-turn  In  one  direction  to  the 
other  in  less  than  1  second.  This  yields  the  constraint 

UtI  -  2  I+tIl- 

Figure  3*2  shows  acceleration  limit,  turn-rate  limit,  and 
the  bank  angle  for  a  coordinated  turn,  as  functions  of  V/Vg 
where  Vg  Is  the  stall  speed  for  a  fighter  in  combat  configura¬ 
tion  and  Vg  =  200  ft/sec.  is  assumed  as  a  nominal  value.  The 
relationships  for  a  coordinated  turn  are 


Acceleration  Limit  -  2.iJ5  Vg 


g' s 

Bank  Angle 
Roll  Rate 


6.0 


,  V-  >  2.^5  V, 


UtIl  “  ^  (l/N^). 


Turn-Rate  Limit 
Turn-Acceleration  Limit 

Iv^Il  -  1-0  g 


=  2  -  2  cos’^  > 

I^T^L  '  Tv/Vg ) 


( 


v/vf  I^tI 


yIl  =  2|YtIl 


Airspeed  change; 
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Equations  3.17  and  Figure  3.2  thus  represent  the 
statement  the  constraints  which  form  part  of  the 
mathematical  model.  In  an  air-air  combat  situation,  it 
will  usually  be  the  case  that  the  target  will  be  flying 
at  a  relatively  high  speed,  so  that  the  left  of  Fig.  3.2, 
representing  the  condition  -  2.^5  V^,  may  be  disregarded, 
and  only  a  relatively  narrow  range  of  V^p/Vg  can  be  realisti¬ 
cally  expected. 

TT  R/S  <t>° 


Figure  3.2.  Bank-Angle  and  the  Limits  of  Target  Turn-Rate  and  Acceleration  versus 

V/Vj;  (Vs  =  Vstall  =  200  VS) 
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A  NONLINEAR  ESTIMATOR 

A  derivation  of  a  nonlinear  range/range-rate  estimation 
filter  Is  presented  below. 

This  filter  has  a  format  such  that  it  is  very  easy  to 
impose  and  take  advantage  of  the  physical  constraints  which 
limit  the  maneuverability  of  an  airborne  target.  It  there¬ 
fore  could  be  appropriately  used  alone  if  initialized  or 
preset  at  launch  with  range  and  range-rate;  thereafter  it 
will  correctly  keep  track  of  those  variables.  Alternately 
it  may  be  used  as  a  nonlinear  filter  in  conjunction  with 
other  ranging  r.cthods,  v.'hether  direct  or  indirect ,  The 
analysis  below  assumes  that  this  estimator  w.ill  be  initialized 
or  preset  at  launch.  The  error  analysis  and  performance  of 
Sections  4.2  and  4  3  show  the  effects  of  errors  in  Iriitializing 
as  well  as  the  effects  of  various  other  error  sources. 

This  nonlinear  filter  has  as  its  physical  origin,  and 
therefore  motivation,  the  hypotheses  that  (a)  the  target  cannot 
easily  change  airspeed,  and  (b)  range-acceleration  is  not 
trivial  and  must  be  estimated. 


4.1  Theoretical  Basis 


The  equations  which  present  the  components  of  velocity 
parallel  and  perpendicular  to  the  LCS  are,  respectively,  from 

(3.15)  and  (3.16), 


r 


V 

T 

r 


ro 


(4.1) 

(4.2) 
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where  the  subscripts  (T)  and  (M)  refer  to  target  and 

missile,  respectively,  while  the  subscripts  (r)  and  (o) 

refer  to  the  components  parallel  to  the  LOS  '’’pendicu- 

lar  thereto.  Thus,  for  example,  the  targ^.  velocity 

component  perpendicular  to  the  LOS  is  . 

0 

Similarly,  the  equations  which  present  the  components  of 
acceleration  parallel  and  perpendicular  to  the  LOS  are, 
respectively,  from  C3.7)  and  (3.10), 

.  2 

r  -  ro  -  a„  +  ajy[  =  0  (I4.3) 

r  r 


ro  +  2fo  -  +  aj,.  *=  0,  (^1.4) 

■*■0  o 

wher..  vhe  components  use  the  same  subscript  code  described 

above,  so  the  target  acceleration  across  the  LOS  is  a^p  . 

Consider  Eqn.  (^.3);  we  can  measure  o  direct-/  from  the 

2 

seeker  and  can  therefore  compute  ro  ,  as  initial  values  of 

r  and  r  are  assumed  given.  We  can  also  directly  sense  the 

missile  acceleration  component  a.,  .  If  we  could  also  estimate  a„ 

we  could  compute  r  and  then  by  Integration  continue  to  estimate 

f  and  r.  This  v/oulcl  enable  us  to  close  a  computation  loop 

and  keep  the  estimation  process  going.  It  v.'lll  be  shown  that 

we  can  estimate  a^  by  using  the  assumption  that  target 

^  r 

airspeed  is  quasi-const  .nt  during  the  brief  interval  of  an 
air-air  encounter.  This  is  significant,  as  linear  estimators 


uniformly  fail  to  operate  successfully  whenever 
negligible,  as  will  be  seen  in  Section  . 


is 


not 
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Let  us  solve  Eqns.  and  (M.2)  for  the  target 

velocity  components  V™  and  V_  .  We  know  that  the  total 

■‘r 

velocity,  or  airspeed,  ,  of  the  target  can  be  de¬ 
termined  from  these  components  by  .>e  relationship 


due  to  the  definitions  cos  and  sin  C,j,. 

r  c 


An  air-air  combat  target  cannot  easily  change  its 
airspeed  V^.  Since  is  nearly  constant,  Is  nearly 
zero,  so  that  differentiation  of  (4.5)  yields 


(4.5) 


■t  Vfp 

r  r  0 


:  0 


(4.6) 


where  a^  and  a,p  are  the  components  of  target  acceleration 
r  o 

parallel  and  perpendicular  to  the  LOS,  The  considerable 
algebraic  effort  required  to  show  that  tnls  differentiation 
is  valid  is  presented  In  Appendix  A. 


We  may  solve  (4.6)  for  a^  ,  so  that  we  may  estimate 
_  _  0  0 


(4.7) 


We  now  substitute  for  a^  >  V„  , and  from  (4.1),  (4.2) 

■^o  '^0  "^r 

and  (4.4)  so  that  (4.7)  becomes 


(ro  +  (ro  +  2^0  +  a^^^ 


(4.8) 
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All  quantities  on  the  right  of  (^<.8)  may  be  measured  or 

At 

estimated  in  the  missile,  so  that  we  can  compute  a,p  .  But 
we  also  ha'i-e,  from  Section  3,  Eqn.  (3-8) 


sin  Krj,  +  Yrpcos  €,p 
o 


and  if  V„  ~  0,  then 


(^.9) 


o  r 

We  may  therefore  compute  a,p  in  a  different  way,  as 

r 

/s 

A  .  Ym  COS  5m  ,  ^  A 

V- 


(^.10) 


where  is  the  estimate  of  y^  »  target  turn-rate. 


Eqns.  (^.8)  and  (^.10)  present  different  forms  of  the 
same  result;  either  may  be  used  as  appropriate  to  the  availa¬ 
ble  physical  constraints.  In  particular,  the  relationship 


(^.11) 


may  be  used  to  estimate  Yip- 

If  we  substitute  (^J.ll)  into  (^.7)  we  have 


ap  »  -Vp  Yp  ^ 

r  a 

This  is  an  obvious  result:  Vp  is  the  target  velocity 

across  the  LOS,  and  Yp  is  its  turn-rate,  therefore  the 
product  is  the  acceleration  of  the  target  in  the  range 


(^.12) 
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direction.  Its  magnitude  is  at  a  maximum  when  the  target 
is  flying  perpendicular  to  the  LOS  and  vanishes  in  a  head-on 
or  tail-chase  combat  configuration. 

Equation  (^.11)  is  a  very  convenient  result,  for  it 
enables  implementation  of  the  constraints,  described  in 
Section  3,  on  the  target  maneuver  capability; 


(1)  The  estimates  a„  and  a,p  may  be 

■^r  a 

limited  to  6g  magnitude. 

(2)  Since  may  be  estimated  as  Y,p  , 


it  is  possible  to  lln:lt  the  quotient  to  a 
realistDC  turn-rate  v.'hlch  may  be  either  a 
nominal  constant  or  a  function  of  the  target's 

A 

estimated  airspeed  V^. 

(3)  Further,  If  Y.p  is  formed  in  a  rate-limited 
net,  it  is  possible  to  Impose  a  physically 
motivated  constraint  on  to  represent  the 

finite  time  required  by  the  target  to  change  its 
turn-rate,  e.g*,  a  right  turn  to  a  left  turn. 
These  limit-properties  are  of  great  value,  for 

'i 

they  make  it  possible  to  estimate  k'lthout 

these  limits,  the  division  in  (^.11)  would  fail 
due  to  a  "division  by  zero"  whenever  the  target 
velocity  vector  is  nearly  perpendicular  to  the 
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LOS.  This  condition  can  occur  in  a  dogfight 

A 

situation.  Further,  the  limit  on  Vtj.  noted 
above  in  (3)  enables  exclusion  of  noise  which 
implies  impossible  roll-rates,  Just  as  the 
limits  (1  &  2)  preclude  undue  sensitivity  to 
noise  implying  impossible  turn-rates  and 
accelerations . 

The  result  of  the  analysis  above  is  the 

estimator 


Fig.  iJ.l  shows  an  estimator  configuration,  which  re¬ 
quires  6  multiplications  for  its  instrumentation  in  the  given 
format.  The  locations  and  character  of  the  limits  which  impose 
the  physical  constraints  are  shown.  The  adaptive  algorithm 
which  determines  Y,p  is  shown  as  part  of  the  overall  diagram. 

It  is  assumed  that  missile  acceleration  components  are  sensed 
with  respect  to  the  seeker  centerline;  it  is  also  assumed  that 
the  milssile  velocity  components  with  respect  to  the  seeker 
centerline  are  available.  This  net  can  be  formed  in  a  variety 
of  equivalent  ways;  the  configuration  shown  in  Fig.  4.1  is 
representative  but  not  unique. 
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A  brief  discussion  of  Pig.  *1.1  follows.  The  block 
at  the  extreme  upper  left  shows  the  process  of  inserting 
the  initial  values  of  range  and  range-rate  at  launch. 

The  three  blocks  immediately  below  show  the  various 
signals  used  in  this  estimator,  gathered  by  their  several 
sources..  Block  A  shows  the  two  integrations  relating 

A  ^  ^  A 

5^  to  r  and  r.  Block  B  shows  the  formation  of  a^  from 

o 

the  relationship 


A  A  A  ^ 

a.^  *  rJJ  +  2fo+aj^j  . 

■*^0  o 


(il.lU) 


Blocks  C  and  D  show  the  formation  of  V_  and  V„  , 

■^r 

respectively,  from  the  relationships 


r  +  V 


(^♦.15) 


ro 


+  V, 


M. 


(i>.l6) 


Block  E  shows  the  formation  of  f  as 


(i».17) 


We  have  stated  that  the  term  a^  is  estimated  as 


^T  "‘^T  ^T  (4.1 8 
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and  chat  this  may  be  limited  to  6g,  or  some  similar 
limit  representative  of  the  abilities  of  the  target 
airplane  and  pilot.  This  limited  product  is  contained 
within  Block  P. 

We  now  consider  estimation  of  Y.p  »  noting  that  if 
V.J,  =  0,  then 


We  do  not  have  the  true  values 
their  estimates.  The  procedure 
form  an  error 


of  arp  and  ,  but  onl 

for  estimating  y  is  to 


y 


e 


Yt^T 

r 


(^.19) 


as  shown  in  Block  G.  A  least  mean  magnitude  algorithm 
is 

A 

Yrp  =  kV„  sgn(e)  (i4. 20) 


where  sgn(e)  is  defined  as 


sgnCe ) 


e  >  0 
e  <  0 


i.e.,  "the  sign  of  c".  One  simple  way  to  limit  y.j 
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to  the  values  permitted  by  the  constraint  relationship  is 

to  replace  by  sgn(V_  )  ,  and  to  choose  k  eaual  to 
■^r 

the  constraint  limit.  It  is  essential  that  the  algorithm 
contain  at  least  the  sign  of  and  the  sign  of  e  if  it 

ys.  X*  /^. 

is  to  be  stable  for  Y,p  >  0  and  also  for  Vrp  <  0.  With  these 
comments,  the  algorithm  is  now 


=  k  sgn(V^  )  sgn{£ )  . 
r 


(4.21) 


The  least  mean  magnitude  algorithm  is  preferred  to  the 
least  mean  square  as  it  yields  more  rapid  solution  for 
this  class  of  problem,  and  is  less  complex  to  instrument. 
In  some  extreme  trajectories  the  target  is  flying  nearly 
perpendicular  to  the  LOS  and  V„  niay  be  quite  small. 

■^r  j. 

To  preclude  instability  of  the  Ym  algorithm  it  is  useful 

/N  A 

to  disable  integration  of  Vrp  when  |V.p  \  is  less  than  some 
minimum  value,  such  as  60'/sec.  Assume  a  target  airspeed 
of  *  lOOO'/sec.  This  limit  is  6"?  of  the  airspeed, 
implying  that  the  angle  between  the  LOS  and  tlie  target 
velocity  vector  is  in  the  range  86®  <  l6n-,t  94°. 

With  this  restrlctloi, ,  algorithm  (4.21)  becomes 


k  sgn(V^  )  3gn(c)  ,  | |  >  V, 


,  otherwise 


^  min  * 


(4.22) 


An  alternate  form  of  (4.22)  which  contains  the  entire 


algorithm  in  one  expression  is 
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A 


k  sgn(V^  )  sgn(e) 


(4.22a) 


We  found  that  k  =  0.4  and  V_  =  60  ft/sec  are  satis- 
factory  values  for  the  missile  and  target  which 

we  assumed.  Eqn.  (4.22)  is  Instrumented  in  Block  H. 

A  yv  A 

The  integration  of  Vrp  to  form  y.j.  and  the  limitation  of 
to  the  constrained  range  appear  in  Block  K.  The  expression 
in  the  square  brackets  [  ]  in  (4.22a)  introduces  the  dead- 
space  in  the  algorithm  and  is  visible  in  Block  H.  The 
nonlinearity  in  the  feedback  path  in  Block  K  is  an  analog 
representation  of  a  saturation  limiter  which  Imposes  the 
constraint  Y-p  on  turn-rate  estimate  Ym.  The  net  of 

Li  ^ 

Fig.  4.1  can  be  mechanized  in  analog  or  digital  formats 
with  equal  ease.  The  choice  of  an  analog  format  for  this 
Figure  is  arbitrary. 
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Figure  4. 1 .  Range  and  Range-Rate  Estimator  Configuration 
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^■2  Error  Analysis 

The  true  value  of  the  range  acceleration  is 

.  2 

r  =  ro  -  a  a  . 

r  r 

and  we  have  shown  that 


^T 

cos 

Cip  "  ^ fjiY ip 

sin 

=  (v^v^-v^  a^  )/V„ 

r 

■^o  ■^o  "^r 

where  = 

■^r 

^T 

cos 

^T  " 

cos 

C 

►a 

1 

II 

sin 

Crp  =  ro  + 

% 

sin  £„ 

o 


and 


~  ’’o+2ra  +  . 


Disregarding  the  effects  of  target  velocity-change  V^, 


the  algorithm  estimates  as 

■^r 


(ro+^M  sin^M) 
r+Vj^  cos 


(ro  +  2  ro  + 


imposing  limits  on  j  a,j,  |,|Y,pl»  and  |V^1  in  the  process. 

r 

We  consider  in  this  section  the  several  sources  of 
error  and  their  significance.  The  principal  sources  of 
error  are  discussed  below. 


(^.23) 


(1)  V^.  The  target  may  change  its  airspeed,  although 
relatively  slowly. 

A  A 

(2)  r^,  r^.  Range  and  range-rate  may  be  initialized 


with  errors. 
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.23) 


(3)  Vj^.  The  missile  airspeed  is  required.  It  may  be  sensed 
or  may  be  provided  adaptively,  but  with  an  error  in 
either  case. 

(^)  ^j,..  The  missile  angle  of  attack  must  be  estimated 
to  form  *  -  a  -  X  in  order  to  resolve  the  missile 
velocity  into  the  components  Vj^  sin  and  Vj^  cos 
Error."^  enter  as  incorrect  estimates  of  A  and  a,  as 
discussed  in  the  simulation  results.  The  angle  X  is 
the  angle  from  the  m.lsslle  centerline  to  the  seeker 
centerline . 

(5)  o.  The  LOS  rate  may  have  random  noise  errors  plus 
a  bias. 

(6)  I  Sjvj  #  The  accelerometers  sensing  missile  accelera- 
r  0 

tion  parallel  and  perpendicular,  respectively,  to  the 
LOS  may  have  bias  errors  and  random  noise. 


We  define  x=r-r,  x=r-r,x=r-r, so  that  x  is  the  range  error. 


and  V.  =  V.  +  N.  ,  . 

ill  v.'hers 


^1  = 


=  0;  Nj  =  shows  Vrr,  /  0 

V  =  V 
3  M 


=  0 

Vg  =  a 

"7  = 
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Thus,  represents  the  deviation  from  the  actual  value 
of  V^,  i.e.,  the  additive  noise.  The  error  x(t)  is  then 
governed  by  the  equation 

X  +  a,(t)x  +  ap(t)x  =  J  b.(t)N,  (4.2^1) 

where  the  coefficients  a^,bj^  are  defined  ?s 
=  2a  tan  -  Y.p  tan  ^ 

^2  =  +  YipCi  +  a  tan  C,p 

bl  = 

=  sec 

^3  "  “^T  ^T 

b^  =  -^Yrp  -2r  tan  Cj.  +2ro 

bg  =  -r  tan  C-j, 
by  =  -tan  . 

The  coefficients  of  x  and  x  are  functions  of  time.  The 

trajectory  of  the  errors  of  rar.ge  estimate,  x,  and  range-rate 

estimate,  x,  may  be  found  by  solving  this  linear  differential 

equation  v/ith  time-varying  coefficients,  subject  to  the  Initial 

conditions  x  and  x  ,  which  are  the  assumed  initialization 
o  o’ 

errors  of  range  and  range  rate.  As  (t.2^)  is  linear,  super¬ 
position  is  valid  so  that  it  may  be  solved  for  the  Individual 

terms.  The  coefficients  a,  and  b^  are  tabulated  in  Appendix  B 
for  a  number  of  typical  trajectories. 
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4.3  Performance 

The  alr-alr  missile  homing  situation  was  programmed 
on  the  digital  computer.  The  simulation  used  parameters 
from  a  typical  short-range  highly  maneuverable  alr-alr 
missile.  The  guidance  method  used  for  the  simulation  was 
the  conventional  proportional  guidance; 

"c  =  5  a, 

where  a  is  the  LOS  angular  inertial  rate  and  s,^  is  the 
acceleration  command  across  the  LOS.  The  transfer  function 
describing  the  missile  characteristics  was 


a(s)  _ 


-53s^-t-106. 6s•^27Q0Q 
s^+72s^+21605+27000 


where  a  is  the  cross-body  acceleration.  This  corresponds  to 
a  simplified  and  linearized  second  order  model  for  the  missile 
airframe  combined  with  a  first  order  model  for  actuator  lag. 

The  estimates  of  range  and  range-rate  were  thus  not  used  in 
the  guidance  loop  as  they  might  be  in  an  actual  missile  system. 
The  missile  was  started  with  a  specified  velocity  with  the 
airframe  state  variables  at  rest  at  the  beginning  of  the  tra¬ 
jectory.  The  geometry  of  the  combat  situation  is  shown  in 
Fig.  4,2. 

A  discretization  seep  size  of  .005  seconds  was  found  to 
result  in  negligible  error  in  the  numerical  output  of  the 
digital  program  as  a  substitute  for  the  actual  continuous 
airframe  being  analyzed. 
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Several  different  geometrical  situations  were  simulated, 
corresponding  to  different  angular  orientations  and  bank  angle 
of  the  target.  Missile  velocity  was  Mach  2.0  or  21^^  ft/sec 
at  10,000  ft.  altitude,  and  target  velocity  1,000  ft/sec.  The 
missile  airspeed  was  assumed  constant;  this  does  not  restrict 
the  generality  of  the  results. 

The  simulations  which  follow  are  arranged  in  ten  groups, 
each  corresponding  to  a  different  starting  situation.  The 
angle  (t.p  is  defined  as  the  bank  angle  of  the  target.  For 
example,  80°  represents  a  hard  left  turn  of  approximately 
^.7  g's  lateral  acceleration.  The  angle  Yrp  represents  the 
target's  initial  heading,  measured  positive  counter-clockwise 
from  a  reference  heading  parallel  to  the  attacking  missile's 
initial  centerline,  which  is  the  x-axls. 

Most  of  the  ten  groupings  of  plots  have  six  separate 
graphs,  arranged  in  pairs  on  facing  pages. 

The  first  graph  is  a  view  looking  down  at  the  two  dimen¬ 
sional  chase.  In  each  case  the  target  is  intercepted  as  a 
result  of  the  homing  navigation.  This  also  shows  initial 
physical  orientation  of  the  target.  Connecting  the  missile 
and  target  at  corresponding  successive  instants  of  time  are 
dashed  lines  representing  the  position  of  the  line-of-sight 
(LOC)  as  the  chase  evolves. 

The  second  graph  shows  the  actual  values  of  range, 
range-rate,  and  range-acceleration  as  functions  of  time. 
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fror.  t  =  0  to  t  =  time-of-lntercept .  These  variables 
represent  the  true  values  of  the  parameters  which  are 
estimated. 

Graphs  3  and  ^  in  the  group  plot  the  errors  In  the 
range  estimate  due  to  various  error  sources.  The  plot 
lettered  (a)  has  no  error  sources  and  represents  the  good 
estimates  of  range  and  range-rate  obtainable  with  this 
method.  The  reason  (a)  is  not  identically  zero  arises 
from,  the  Inaccuracies  of  the  algorithn.  ana  computation;  it 
is  intended  as  a  control  in  analyzing  the  other  graphs. 

The  other  plots  show  the  effects  of  various  errors.  For 
example,  in  plot  (b)  range  is  initialized  incorrectly  at 
99%  of  its  actual  value.  A  listing  of  the  sources  of  error 
is  as  follows: 

(a)  No  errors,  except  those  inherent  to  the  algorithm 
Itself. 

(b)  Range  initialized  with  1%  error. 

(c)  Range-rate  initialized  with  1%  error. 

(d)  The  signal  o  is  erroneously  biased  by  the 
amount  0.5  deg/sec. 

(e)  The  signal  G,  the  lateral  accelerometer 

p 

output,  is  biased  by  3  ft/sec  . 

(f)  The  signal  H,  the  longitudinal  accelerometer 

2 

output,  is  biased  by  3  ft/sec  . 

(g)  Tl;e  seeker  angle  relative  to  the  missile 
centerline , X ,  is  biased  by  0.1  degrees. 

(h)  The  error  in  measuring  missile  velocity 
is  taken  as  1%  or  21.^1^  ft/sec. 
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(i)  The  error  in  estimating  angle  of  attack, a, 
was  taken  to  be  a  bias  of  1  degree. 

(J)  In  this  case  the  target  velocity  is  increased 
with  an  acceleration  of  Ig ,  thus  illustrating 
the  effect  of  violating  the  assumption  V^=0. 

In  all  cases,  only  the  particular  source  of  error  mentioned 
was  included,  i.e.,  the  other  error  sources  vjere  removed. 

Graphs  3  shows  the  errors  from  sources  (a)  through  (e) 

above.  Graph  ^  shov;s  error  sources  (f)  through  (j). 

In  the  case  of  Figure  ^.6,  an  additional  source  of  error 
is  analyzed.  Figure  shows  the  effect  of  additive  noise 

in  the  LOS  rate  as  used  by  the  kinematic  ranging  algorithm. 

The  corrupted  LOS  rate  was  also  fed  to  the  guidance  loop. 

Graphs  5  and  6  of  each  group  show  the  corresponding 

errors  in  range-rate  instead  of  the  error  in  range.  Separation 
of  cases  (a)  through  (J )  is  the  same  as  above  in  graphs  3  and 

4. 

Figure  4.12  is  also  an  exception  to  the  format  specified. 
This  figure  displays  the  geometry  when  the  target  changes  its 
maneuver  during  the  flight  of  the  missile.  The  target  starts 
with  a  turn  to  the  left,  changing  to  a  turn  to  the  opposite 
direction  of  equal  magnitude. 

Figure  4.12.2  shows  both  the  error  in  range  and  range- 
rate  from  the  estimation  algorithm.  This  error  represents  the 

numerical  accuracy  of  the  algorithm,  and  is  not  due  to  other 
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FIGURE  4.3 

Graph  1,  the  upper  graph  to  the  left,  shows  the  plan  view  of 
the  geometry  of  the  missile  chase.  The  target  has  an  initial  range 
of  slightly  more  than  10,000  ft.,  and  initially  is  traveling  across 
the  LOS.  The  target  airplane  undergoes  a  5.7g  constant  right  turn 
which  is  indicated  by  the  bank  angle  =  -80°. 

Graph  2,  bottom  left,  plots  actual  range,  range-rate,  and  range-  ! 
acceleration.  One  may  note  the  decrease  in  magnitude  of  range-rate 
as  the  chase  changes  fromi  a  side  attack  to  a  tail  chase. 

Graphs  3  and  4,  top  and  bottom  immediate  right  on  the  facing 
foldout  page,  show  the  errors  in  the  range  estimate  due  to  various 
error  sources  for  the  geometry  described  above. 

Curve  (a)  represents  no  errors  except  those  due  to  the  comiputatlon 
process  and  shows  the  ability  to  track  the  target  when  perfect  signals 
are  used  in  the  algorithm.  Maximum  error  for  range  is  approximately 
5  ft.  Curves  (b)  thru  (j),  on  Graphs  3  and  4  represent  the  effects 
of  the  error  sources  specified  in  Section  4.3. 

Graphs  5  and  6,  on  the  far  right  of  the  facing  page  display  the 
errors  in  range-rate  from  the  error  sources.  Curves  (a)  thru  (J) 
again  correspond  to  the  errors  mentioned  in  Section  4.3.  Curve  (a) 
shows  a  maxim.um  error  of  3  ft/sec. 

Curves  (J )  on  Graphs  4  and  6  represent  the  largest  errors,  and 
appear  reasonable  since  the  basic  assumption  pf  the  algorithm  (V^=0) 
is  violated.  Curve  (J )  thus  shows  the  error  due  to  target  airspeed 
changes , reading  a  maximum  range-error  of  350  feet. 

For  convenience,  the  labels  a...J  in  the  curves  imply: 


no  error  sources 
\%  error  in  r  (0) 
1%  error  in  r  (0) 
0, 5°  bias  in  o 
0.1  g  bias  in  a„ 


C.l  g  bias  in  a^^ 

0.1  bias  in  seeker  gimbal  angle 
20  ft/s  bias  in 
1°  error  in  a 
V^.  =  1  E 
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FIG.  4  .4  .1  OYNfiMlC  RANGING  SIMULATION  WITH 
TARGET  INITIAL  CONDITIONS.  ♦  =  90  .0=0  .XT  =10000. 
YT  =  1000.  PLOTS  OP  VARIABLES  DEFINING  GEOMETRY 
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FIGURE  11.4 


This  grouping  of  curves  shows  a  non-maneuvering  target, 
with  initial  heading  90°  relative  to  the  attacking  missile's 
centerline.  In  this  situation,  the  LOS  rate  approaches  zero 
after  the  missile  response  to  the  guidance.  As  a  result,  one 
can  say  in  general  that  errors  in  the  estimates  will  tend  to 
grow  lar^e  after  LOS  rate  gets  small.  This  is  somewhat 
apparent  to  the  curves  on  the  right.  It  is  also  noted  that 
curves  (a)  are  not  small  compared  to  the  others,  as  was 
observed  in  Figure  3-  This  is  explained  by  the  same  reason, 
namely,  that  when  neither  the  target  nor  missile  maneuver, 
the  estimation  of  range  and  range  rate  is  difficult  and  more 
sensitive  to  computation  error.  The  algorithm  exhibits 
moderate  sensitivity  to  initial  range  error  and  to  accelerometer 
biasses . 
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FIG.  4  .H  .2  dynamic  RANGING  SIMULATION  WITM 
TARGET  INITIAL  CONDITIONS.  t=90  ,  ♦=  0  .  XT  =  10000. 
YT  =  1000.  ERRORS  IN  RANGE  ESTIMATE  FROM  VARIOUS 

SOURCES 
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FIG.  4  , 4  .3  DTNRMIC  flflNGING  SIMULATION  WITH 
TARGET  INITIAL  C0N0ITI(3NS.  ♦  =  90  ,  ♦=  0  .XT  =  10000  . 
YT  =  1000.  ERRORS  IN  RANGE  RATE  ESTIMATE  FROM 
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FIG.  4  .4  .1  DYNAMIC  RANGING  SIMULATION  WITH 
TARGET  INITIAL  CONDITIONS.  t=90  .  ♦=  D  ,  XT  =  -0000 
YT  =1000.  PLOTS  OF  VARIABLES  DEFINING  GEOMETRY 
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FIGURC  ‘4.5 


Figure  ^4 . 5  shows  the  target  with  an  initial  heading  of  90*^ 
relative  to  the  missile  centerline,  and  a  bank  angle  of  8o'^, 
implying  a  hard  turn  to  t'ne  left.  With  this  configuration,  the 
aspect  angle  5rj-  is  close  to  90°  at  all  times,  in  fact  it  passes 
through  90°.  The  error  analysis  equation  predicts  that  the 
estimates  would  be  extremely  sensitive  in  this  case.  However, 
the  nonlinearity  in  sensing  the  sign  of  ,  which  is  not 
considered  in  the  error  analysis  equation,  prevents  the  estimates 
from  becomiing  sensitive,  and  In  fact  the  errors  from  all  sources 
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Without  the  dead  space  in  Block  H  of  Fig.  4.1,  the  estimator  would  be 
unstable,  see  Fig.  4.5.4  below. 


IflBGfl  INITIAL  CONDITIONS.  t=90  .  ♦=0O.  xT  =  10000 

TT  =  1000.  ERRORS  IN  RANGE  ESTIMATE  &  RANGE  RATE  ESTIMATE 

FROM  VARIOUS  SOURCES 
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FIGURE  4.6 


This  shov/s  the  target  vilth  an  initial  relative  heading 
of  0°,  and  a  bank  angle  of  80°.  The  largest  error  appears 
to  occur  for  case  (j),  target  airspeed  change,  In  agreement 
with  previously  explained  simulations.  The  other  significant 
error  source,  (d) ,  is  LOS  rate  bias.  » 

FIGURE  4.6.4,  below,  shows  the  errors  v;hen  noise  is  added  to  the 
LOS  rate,  instead  of  the  fixed  bias  present  in  curve  (d).  This 
random  noise  wa?  constructed  to  have  zero  mean  and  standard 
deviation  equal  to  5*  0/Range;  thus  the  noise  becom.es  quite  large 
as  range  decreases. 


a.  no  error  sources  f.  0.1  g  bias  in  a^^ 

0  ^ 

b.  1^  error  in  r  (0)  g.  0.1  bias  in  seeker  gimbal  angle 
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FIG.  li  .6  .3  DYNAMIC  RANGING  SIMULATION  WITH 
TARGET  INITIAL  CONDITIONS.  t=0  .  ♦=80.  XT  =10000  . 

YT  =  1000.  ERRORS  IN  RANGE  RATE  ESTIMATE  FROM 
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FIGURE  ^.7 


This  geometry  displays  a  head-on  situation  where  the 
target  and  missile  are  approaching  each  other  with  maxim.um 
closure  rate.  The  target  Is  turning  left.  Error  curve  (a) 
remains  small.  Curve  (b)  is  worthy  of  note,  as  it  demonstrates 
that  an  error  in  initial  range  estimate  need  not  lead  to 
divergence  of  the  range-rate  estimate.  Curve  (c)  shows  readily 
that  an  error  in  range-rate  leads  to  an  increasing  error  in 
range.  Curve  (J  )  again  shows  the  largest  errors,  due  to  target 
airspeed  change. 
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FIGURE  <4.8 


The  figure  shows  a  head-on  situation,  with  the  target 
initially  at  slightly  more  than  20,000  feet  range,  flying 
straight  (no  turn  rate).  Range-rate  is  quite  constant,  and 
range  acceleration  is  negligible  and  barely  discernible 
from  the  time  axis.  Errors  in  the  estimates  are  very  small 
except  for  a  fev;  cases.  Curve  (c)  again  shews  that  a  con¬ 
stant  range-rate  error  leads  to  an  increasing  range  error. 
Again,  the  range  error  in  curve  (b)  does  not  lead  to  a  range- 
rate  error.  Curves  (j )  show  the  largest  errors. 
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FIGURE  tj.9 


Thj.s  figure  has  tarp-c-t  headinf  0°,  and  bank  angle  -80*^. 

Thus  pursuit  starts  as  a  tail  chase,  and  ends  with  a  rear-side 
attack.  Target  initial  range  is  sliglitly  more  than  10,000  ft. 

The  estimate  errors  generally  fellow  the  sarrie  pattern  established 

in  [-revious  trajectories,  with  the  largest  error  sources  being 

« 

0 


bias 

(d)  and 

target  airspeed 

change 

(j). 

a. 

no  error 

sources 

A 

f.  0.1  { 

'  bias  in  a,. 

% 

b . 

1%  error 

in  r  (0) 

g.  0.1° 

bias  in  seeker 

c . 

Irt  error 

A 

in  f  (0) 

A 

h.  20  ft/s  bias  in 

d . 

0,5*^  bias 

« 

in  0 

i.  1°  error  in  a 

e . 

0.1  g  bla 

s  in  a,. 

J  •  v^.  = 

1  g 

o 

o 


fIG.  4  .10.1  LnrNflMJC  RANGING  SIMULATION  WITH 
TflRCETT  INITIAL  CONdlTIONS.  +-180.  ^.-r-80,  XT  ^  10000, 
TT  1000.  PLOTS  or  VARIABLES  DEFINING  GEOMETRY 


Bell  Aerospace  Company 


FIGURE  4.10 


This  trajectory  displays  a  head-on  trajectory  with 
initial  range  around  10,000  ft.  The  target  is  executing 
a  hard  right  turn.  The  estimate  errors  are  somev.-hat  re¬ 
lated  to  the  errors  portrayed  in  Figure  4.7,  where  the 
turn  was  in  the  opposite  direction. 
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FIGURE  ^.11 

Figure  ^.11  portrays  a  side  attack.  Initial  range  is 
slightly  more  than  20,000  ft.,  and  the  target  is  not 
maneuvering.  As  is  evident,  the  guidance  system  has 
established  an  intercept  course  at  around  3  seconds,  and 
after  this  point,  missile  maneuvering  is  negligible.  This 
results  in  the  signals  o  and  consequently  a  approaching 
zero,  and  thus  the  main  signals  being  used  by  the  algorithm, 
have  little  significance.  Thus  after  3  seconds  it  can  be 
ar.ticipatcd  that  if  the  algorithm  is  not  set  up  properly  at 
that  point,  then  serious  errors  can  and  will  result.  This 
is  evident  in  the  curves  to  the  right. 
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FIGURE  4.12 

displays  a  sltuatj.on  v.'liere  tU.e  target  pci-forms 

a  rarlvi  cr aiire  In  acceleration  by  changing  its  turii  from  left 

'c  :  .  This  trajectory  is  identical  to  that  of  Fig.  ^.6  for 

Iv  first  three  seconds ;  at  t  =  3,  the  target  bank  angle  was 

cb'.-.nge-d  from  8o'"'  left  turn  to  8C°  right  turn,  requiring  one 

second  to  do  so;  the  bank  angle  is  80°  right  turn  for  t>^.  The 

al.^oriti.:.  v.t.ich  estimates  range  must  therefore  re-adagt  to  acconi- 

r.odater  the  chaiige  in  target  turn  I'ate.  The  bank  angle  was  changed 

in.  a  linear  fashion;  hiOwever  as ,  acceleration  is  related  to  tan  9, 

tr.c  acceleration  of  the  target  changed  in  a  nonlinear  fashion. 

i;.i  -  Is  evident  or.  inspection  of  the  r  trace,  on  the  left,  at 
<  < 

F  -  t  -  seconds.  As  the  acceleration  changed  in  a  way  v;hich 
violates  t!ie  cst.imatcr '  s  assumptions,  an  error  is  formed  while  the 
:  .ai.'.uv<.r  charige  occurs.  but  there  are  many  ways  in  whlc'  an  air'- 
plane  can  change  its  maneuver,  and  any  one  is  as  likely  as  another. 
The  estimator  properties  lie  in  the  mid-range  of  the  likely  types 
cl'  iiarieuvcr  change  methods.  The  curves  on  the  facing  page  are 
tlic-  error.s  in  the  range  and  range-rate  estimates  with  range  error 
o.M  the  upper  grapli,  and  range-rate  error  on  the  lower.  IJo  addi- 
tlotia]  curves  I’er'rosent ing  sources  of  error  v/ere  made. 

Thi.'  ti'ace  a'  tlie  lower  right  shows  the  error  of  range-rate 

(•;' 1 1  s.at'- .  fcticc  that  this  error  Jumps  from  ft/sec  to  6  ft/sec 

at  t:,'/  tii:,e  (.  f  li.i  '  t  at  the  end  of  the  run.  This  change  is  due 

.2 

to  li'iei’l'ect  tracking  of  the  lerm(ro  )  since  0  is  extreinely  large 
at  i::.|  act  v.'hlie  r  is  not  exactly  equal  to  r;  In  tills  Cotmecliun 

.  2 

.  t!,'.-  d  i  1 01.  ol'  the  physical  meaning  ol'  the  term  ( ro  )  on 


1  aip'  j  : 
•/<-;■  d  a  I  ‘ 
j  I ,  c :  ■ ' ;  a .  ■ 


.!cc‘  also  tliat  the  error  of  tracking  the  chatiglng  maneu- 

<"  < 

L  i  -  t  -  s  seconds  causes  the  I’unge-rate  estii'jate  error  to 
froii.  1  ft. /sec  at  t  =  3  to  ^  ft. /sec  at  t  -  l4;  the  Increase 
of  3  I’t./sec  Is  considered  to  be  mlnliral  and  detronsti'ut.';. 


r  I  : .  ,  n  ('  j  '  .v,  r  o  1'  t  h  !  s- 


nonlitieai-  estlinator. 
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LINEAR  ESTIMATORS 

Linear  estimators  may,  for  the  purpose  of  this  dis¬ 
cussion,  be  divided  into  two  classes,  dependin/5  on  whether 
a  gradient  technique  is  used  to  achieve  convergence.  We 
term  "adaptive"  those  methods  which  use  a  gradient  technique 
in  some  form.  A  matrix-inverse  problem  is  always  implicit 
in  the  formation  of  an  "adaptive"  or  gradient-seeking  solu¬ 
tion.  In  the  "matrix-inverse”  technique,  that  inverse 
problem  appears  explicitly  and  must  be  treated  explicitly. 

The  adaptive  methods  lead  to  simpler  instrumentations 
than  the  matrix-inverse  methods  and  were  therefore  considered 
first.  The  adaptive  methods  failed,  whether  with  simple 
single-error  nets,  or  with  multiple-error  nets  or  Kalman 
filters,  principally  due  to  target  acceleration  effects. 

The  mode  of  failure  was  that  the  estim.ated  range-rate  was  as 
much  as  30%  to  50%  in  error  at  the  end  of  the  major  maneuvers 
of  the  missile.  The  study  then  turned  to  matrix-inverse 
methods;  these  failed  to  converge  for  some  encounter  geometries 
and  converged  to  erroritous  solutions  for  other  geometries.  The 

modes  of  faiJui-e  ■.-.■L-xt.-  lai'ge  <.;:'’rors  of  range-rate  cstl;..'ite  or 
a  probl'-ri'  of  1 1  l-co.'idi  I 'l oiling  v;hlch  pu-ecluded  convei'gence  .  We 
tti'  T''- I’CjI'o  tentalively  conc]ud.;o  that  fundamental  questions  of 
t,,athe!r:at  1  e 0  j  obsorvahi  lily  were  Involved  and  then  verified  that 

I:, VI  otbesis  . 
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This  section  presents  a  brief  outline  of  the  adaptive 
methods,  and  of  the  matrix-inverse  methods,  with  comments 
on  the  physical  reasons  for  their  deficiencies,  followed  by 
a  discussion  on  mathematical  observability,  as  applied  to 
this  situation.  The  section  concludes  with  some  comments 
on  the  possibilities  available  by  combining  the  kinematic 
ranging  linear  estimator  with  other  linear  devices,  in 
particular  with  a  stadlmetric  ranging  method. 

3.1  Adaptive  Estimators 

Adaptive  Identifiers,  as  a  class,  share  the  use  of  the 

gradient  of  the  index  of  performance  as  the  guide  to  im.prov- 

Ing  the  parameter  estimates  from  one  instant  to  the  next. 

The  tv;o  principal  well-known  methods  of  forming  param.eter 

2 

identifying  nets  are  knoivn  as  the  response-error  and  equa- 

1  ;  3  >  5 

tlon  error  methods.  The  response-error  method  is  known  to 
be  relatively  slow,  and  relatively  insensitive  to  noise,  as 
compared  to  the  equation-error  method.  The  air-air  Intercept 
problem  clearly  requires  rapid  solution;  efforts  were  therefore 
concentrated  on  th.e  equation-error  method. 

A  brief  analytical  outline  of  the  method  is  followed  by 
discussion  of  the  representatives  of  tlie  various  adaptive 
systems . 

Assume  that  the  outi)Ut  of  the  plant  whose  para'  eters  are 
to  be  Identii'led  is  designated  y,  and  that  a  set  ol'  sig.ials 
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consisting  of  various  functions  of  the  plant's  input  and 
output,  is  available.  If  the  true  values  of  the  unknown 
parameters  are  represented  by  the  vector  x  then  the  plant 
output  takes  the  form 


Given  the  same  set  of  signals,  v,  and  a  vector  of  estimates 
of  the  parameter  values,  x,  we  may  form  an  estimate  of  the 
plant  output 

y  =  v*^  X .  (5.2) 

Then  we  may  form  an  error  as  y  -  y  =  y  -  v  x  -  e  ,  (5.3) 

and  the  error  w-11  become  identically  zero  if  and  only 

« 

if  X  =  X  .  The  common  and  simplest  gradient  algorithm 

for  adjusting  the  parameter  estimates,  x,  minimizing  J=e^,  is 

X  =  +k  VC  ( 5 . ^  ) 

which  on  substituting  e  from  (5.3)  yields 

T 

(Is+kvv  )x=kvy. 

The  gain  V  must  be  positive  definite;  it  may  be  a  scalar 
or  a  matrix.  If  the  value  of  k  is  chosen  small  enough  to 
ensure  averaging,  and  if  the  parameter  x  is  constant, 
then  the  steady  state  is  reached  v;hen 


T 

v  v  X  =  V  y 

v;lth  nclutlon  x  =  ( - ^  - 

V  V  V  y .  (5.5) 
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The  covariance  matrix  v  v  is  implicit  in  all  adaptive 
techniques  as  they  use  the  algorithm  (5.^)-  The  convergence 
properties  and  noise  sensitivity  of  the  solution  are  de¬ 
termined  by  the  properties  of  that  matrix. 

Matrix-inverse  methods,  to  be  discussed  in  Section  5-2, 

find  the  solution  vector  x  directly  from  (5.5)-  The  co- 

T 

variance  matrix  v  v  is  obviously  Involved  explicitly. 

If  the  parameters  x  are  known  a  priori  to  be  variable 
or  to  be  related,  the  known  properties  may  advantageously  be 
used  in  a  modified  algorithm.  Assume  that  it  is  known  that 
X*  =  x^  and  x|  =  bt,  or  in  general,  x*  =  Ax*  +  B(t),  then  an 
algorithm  equivalent  to  (5.^)  is 


X  =  Ax  +  B(t)  +  k  V  e, 


(5.6) 


and  the  term  k  v  e  is  required  only  to  correct  for  the 
unknown  errors  in  the  parameter  estimates. 

Cur  initial  approach  was  to  use  the  identity 


ro  +  ro  =  y  -  y 


(5.7) 


as  a  starting  point  for  an  adaptive  algorithm.  This  equa¬ 
tion  results  from  the  Inertial  space  geometry 

T 
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which  shows 


r  sin  0  =  -  yj^ 


For  small  o. 


this  equation  has  the  approximate  form 
r  o  =  y^  -  yj^ 


with  first  derivative  (5.7).  Taking  current  estimates  for 
range  and  range  rate,  which  in  general  are  different  from 
the  actual  values,  the  equation  error  is 


XiVi+X2V2+X3V3+yj^  =  e^ 


(5.8) 


where  x^=r^,X2*  r^,X3=  y^,  ~  ~  x^=X2  > 


or  X  =  AX,  A  =  /O  1  0\ 

0  0  0  , 

'0  0  of 

•  « 

where  o,  o  and  y^^  are  known.  The  index  of  error  is  defined 
as 


J  = 


^  T 

The  algorithm  for  updating  the  estimate  x  =  (r,r^y^)  is 
X  =  Ax  -  k  V  e^,  which  is  then  readily  Implemented  for  a 
suitable  choice  of  k  >  0.  Using  this  method  assumes  that  the 

/S  /N 

variables  r  and  y^  are  constant.  But  in  an  actual  dogfight 
situation,  both  target  and  missile  are  accelerating,  and  the 
assumption  is  too  optimistic.  This  approach  therefore  failed 
and  was  abandoned. 
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A  second  approach,  consistent  with  the  mathematical 
model  of  Section  3»  was  taken.  In  this  development,  the 
frame  of  reference  is  polar  coordinates  centered  at  the 
missile,  with  the  instantaneous  radius  vector  parallel  to 
the  line  of  sight.  The  error  equation  which  holds  in  this 
case  is 

/S  yv 

ro  +  2  r  0  -  a^,  +  a^,  =  e^ , 

o  "o 

and  the  error  index  is  then  chosen  as  J  =  e^  • 

^ 

Again,  let  =  r  ,  X2  =  r,  ~  ^1  ~ 

v^  =  -1.  Then  that  error  equation  can  be  put  into  the 

form 

T  , 

V  X  +  a,,  =  e„  . 

% 

The  complete  statement  of  the  components  of  accelera¬ 
tion  is  from.  Section  3, 


(5.9) 


m  i 

V2=o, 

the 

Accelerations 
perpendicular 
to  ohe  LOS 

Accelerations 
parallel 
to  the  LOS 


r  d  +  2  r  a  -  a^  +  =0 

o  a 


r  -  r  o  -  a~  +  a.,  =  0, 

1  I'i 

r  r 


(5.10) 


(5.11) 


The  first  of  these  may  be  used  to  form  the  "observation 
equation",  or,  in  adaptive  terminology,  the  equation  error. 
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But  the  second  expression  provides  a  statement  of  the  value 
of  r  which  has  previously  been  lacking.  Let  us  define 

^  A  /S 

=  r,  X,  =  a^  and  Xj^  =  .  We  can  directly  sense  the 

■^o  ■^r 

LOS  angular  rate  and  acceleration  and  the  unit  gains 

Vt  =  -1,  v^  =  +1,  and  also  the  missile  accelerations  a*,  and 
^  *0 

a„  .  Then  the  equation  error  to  be  minimized  is  again  in 


the  form 


T 

V  X  +  a. 


(5.12) 


and  the  algorithm  which  m.inimizes  J  =  e^  >  is 
x=Ax+B(t)-kve^ 


where 


0 

1 

0 

0 

0 

.  2 

0 

0 

0 

1  ,  B  = 

-^M 

0 

0 

0 

0 

0  ^ 

0 

0 

0 

0/ 

0  1 

(5.13) 


(5.1^) 


The  matrix  A  expresses  the  differential  relationship 
between  range  and  range-rate  on  the  first  line;  it  expresses 
the  definition  of  r  from  (5.11)  on  the  second  line,  and  on 
the  third  and  last  lines  expresses  the  assumption  that  the 
target  accelerations  with  respect  to  the  LOS  are  quasi- 
constant  . 

.  2 

The  term  ro  in  (5.11)  can  easily  be  as  large  as 
IOC  ft/sec%  e.g. ,  when  in  the  initial  geometry  the  target  at 
a  range  of  10,000  feet  has  1,000  ft/sec  velocity  and  is 
flying  perpendicular  to  the  LOS.  This  is  not  negligible. 
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The  kinematic  ranging  estimation  method  expressed  by 
equations  (5.13)  and  (5.1^)  was  examined  for  a  variety  of 
geometrical  conditions. 

The  simple  adaptive  schemes  following  (5.13)  and  (5.1^), 
with  the  adaption-rate  parameter  k  as  a  scalar  constant, 
were  quite  unsatisfactory  as  there  was  no  single  constant 
value  of  k  which  worked  well  for  all  combat  geometries  or 
throughout  any  single  trajectory.  The  difficulty  is  that  a 
large  value  of  k  is  required  to  find  the  correct  values  of 

At  At 

r  and  r  at  the  start  of  the  trajectory,  while  a  low  value 

is  required  to  prevent  undue  sensitivity  to  noise  in  later 

parts  of  the  trajectory.  Further,  the  rate  parameter  k 

should  ideally  be  proportional  to  (77)  so  as  to  accommodate 

o 

both  to  large  and  small  maneuver  geometries.  The  gain  k 
must  thus  vary  with  time  and  also  with  geometry  to  achieve 
satisfactory  results. 


It  was  therefore  necessary  to  replace  the  scalar  adaption- 
rate  parameter  k  by  a  Kalman  Filter  type  matrix,  K,  where  the 
values  of  K  are  time-varying,  according  to  the  well-known® 
rule 

•  —  T  T 

K  =  -K  V  N  V  K 


(5.15 

where  N  is  the  noise  covariance.  Even  so,  while  this 
approach  would  generate  good  estimates  of  r  and  r  for 
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some  trajectories,  it  failed  when  the  trajectory  was  such 
that  a,p  /  0. 


In  the  geometrical  conditions  in  which  this  Kalman- 
fllter  type  algorithm  works  well,  l.e.,  when  the  target 
acceleration  along  the  LOS,  a^  ,  is  negligible,  a  signifi- 
cant  simplification  is  possible.  Consider  (5.11)  and 
"transform-dlfferentlate-lntegrate"®  by  premultiplying  by 


T  s 

>  i=l,2,...,T^/T^,  yielding,  if  r  Z  o, 


TiS 


T^s+ 


/T  s  \  /  T  s 

T  ('■  +  a,,,J  =  r  (.5^  s)  +  r 


a  + 


(T 


»)-  (i^)  *  (5^  ') 


M. 


(5.16) 


For  these  geometries,  it  is  a  reasonably  good  assumption 
that  arp  Is  a_nasi-constant .  It  therefore  follows  that 


T.  s 

= — -y  a„  1  0  Is  valid.  The  result  is  that  the  estimating 

ilS+i 

variable  ^  3  “  ^-p  t)e  omitted,  and  there  is  one  less  un- 

known  parameter.  Equation  (5.16)  may  be  put  into  the  equation 
error  form 


^  ^  +  il  = 
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where 


T^s+1 


2T  s  T  s  .  T 


If  this  transform-integraticn ,  (5-16),  is  instrunented 

for  several  different  values  of  T^,  a  vector  of  errors 

(e^i  ,  e^.,  ...)  n':ay  be  generated.  This  is  one  of  the  several 

1  ^2 

multiple-error  techniques ‘  This  technique  did  net  work  in 

an  adaptive  approach  although  it  was  better  than  single-error 

techniques.  The  solution  usually  involved  an  ill-conditioned 
T 

V  v  covariance  matrix  and  was  therefore  noise-sensitive. 

When  averaging  over  a  sufficiently  large  Interval  was  used, 

the  unknown  acceleration  a.^,  caused  r  to  be  as  much  as  30"  to 

■^r 

5Q%  in  error  in  1  to  2  seconds. 

The  difficulty  is  that  if  the  filter  time -constants  T. 

are  smiall,  then  the  several  vectors  are  almiost 

parallel  resulting  in  an  ill-conditioning  problem  and  a  poorly 

defined  solution,  Cn  the  other  hand,  if  are  large  enough 

to  avoid  ill-conditioning  then  the  invalid  assumiption  that 

r  I  =0  caused  excessive  error  in  the  range-rate  estimate, 

^  r 


63 


Bell  Aerospace  Company 


The  difficulty  common  to  all  the  methods  is  that 
whenever  the  averaging  time  must  exceed  1/^  to  1/2  second, 
to  enable  the  vectors  vjhich  carry  the  information  to  be 
linearly  Independent,  it  is  essential  to  take  into 

account  the  target  acceleration  along  the  LOS.  On  the 
other  hand  an  averaging  time  of  1/2  second  or  less  is 
not  sufficient  as  the  LOS  does  not  rotate  sufficiently  in 

that  interval  to  permit  determining  a  solution,  in  the 

\ 

presence  of  noise,  and  also  to  correct  the  initial  estimates. 
For  these  reasons,  the  various  adaptive  gradient- seeking 
methods  of  parameter  estimation  were  abandoned. 


5.2  Matrix-Inverse  Methods 

7 

Several  schemes  were  derived  from  the  classical  point 
of  view  of  least  mean  square  estimation.  The  general  theory 
of  this  method  can  be  explained  quite  briefly,  as  follows. 
Consider  the  linear  homogeneous  vector  differential  system 

X  =  A(t )x  . 

( 

Assume  the  observed  output  of  the  system  is  the  vector 


y  =  h'  ( t )  X  +  V 

where  H(t)  is  known  and  v  represents  additive  random  noise. 
It  is  desired  to  estimate  x  knowing  only  A(t),  y_(t)  and  the 
.statistics  of  V. 


Sk 
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f 


i 

/ 


Knowing  A(t),  generate  the  state  transition  matrix  » 
for  the  system  so  that  the  state  at  any  time  can  be 
written  as  a  function  of  the  state  at  another  time.  Thus, 


X  (t+c)  =  4  (tf  4,t )  X  (t)  . 


(b.20) 


Many  estimates  can  be  e"'  ,':  acted.  The  one  which  is 
usually  cor"idered  the  most  aeslrable , however ,  has  the 
properly  that  it  is  an  "unbiased"  estimate  of  x  and  further¬ 
more,  it  the  variance  of  the  expected  error  in  the 

estimate.  This  estimate  is  found’  by  computing 


where 


T  m 

fn  ^  m  * 


^(n)  -I  \ 


YlCtp 


\ 


H(t^) 
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(5.21) 


(5.22) 
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The  matrix  Is  thus  a  "continuously  updated  observa¬ 
tion  matrix"  and  the  vector.  is  a  "series  of  output 

observations".  The  expression  (5.21)  is  valid  provided 
the  noise  v  is  Independent  from  computing  Instant  to 
Instant  and  has  zero  mean.  The  construction  of  the 

matrix  T  is  often  a  formidable  task  as  it  requires 
n 

keeping  in  storage  a  large  amount  of  past  data  H(t^)  so 
that  it  can  be  updated  by  the  state  transition  matrix 
4>(ti,tj)  which  is  at  times  also  difficult  to  obtain. 

Apply  this  theory  to,  typically,  the  reduced 
dimensionality  problem  defined  by  (5.16)  and  (5.17); 
assuiae  that  f  Is  constant,  thus  ^2  constant  after 
convergence.  But  x^  »  r  and  is  not  constant.  Now  we 
have,  at  some  specific  Instant,  and  with  only  one  value 
of  1  =  1, 


Vi(ti)xi  +  V2(tj^)x2  +  y^^  '  e^. 

At  the  next  instant  we  have 

But  while  the  value  of  X2  is  Ideally  the  same  at 
these  two  instants,  the  value  of  changes  continuously, 
for  “  r,  and  r  is  the  integral  of  r.  We  may  modify 
(5.24)30  that  x^  is  the  same  for  the  two  Instants.  If  At 


(5.24 


(5.25 
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Is  the  time  between  samples,  we  may  rewrite  (5.2^1)  as 
v^(tT)  (x^(t2)  +  X2  At)  +  V2  (t^)  ^2  ^1  “  ^5  ’ 

which  may  be  rearranged  as 

Vi(Xi)  x^(t2)  +  (v^  (t^)  +  (tj^)  At)  X2  +  =  e^. 

(5.26 

Now  Eqns.  (5.25)  and  (5-26)  represent  two  measurements 
referred  to  the  same  Instant,  and  may  be  solved  Jointly, 

More  generally,  a  larger  sequence  of  such  Instants  may  be 
gathered  In  the  more  general  form  of  (5.23). 

Thin  method  exhibited  several  modes  of  failure. 

Whenever  any  two  of: 

the  missile  velocity  vector 
the  target  velocity  vector 
•-e  LOS 

•  T 

become  col  'ar,  or  0  *  0,  the  matrix  (T  T  )  became  111- 

’  n  n 

conditioned  ^  indeterminate,  yielding  Inaccurate  and  nolse- 

.s  iitive  :  utions.  In  geometries  when  this  mode  of  failure 

dv.  jur,  this  algorithm  has  a  different  mode  of 

failure,  If  the  duration  of  the  stored  memory  exceeds 

1/2  second,  the  neglected  target  acceleration  parallel  to 

the  LOS,  a^  ,  causes  the  estimate  of  r  to  be  to  ^0%  in 

■^r 

error  within  1  to  2  seconds,  as  the  assumption  that  r  =  0  is 
too  Incorrect.  On  the  other  hand,  a  shorter  duration  memory 
results  In  high  noise  sensitivity.  But  these  algorithms  work 
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well  (Range-Rate  errors  of  IJ!)  if  a™  *  0. 

^r 

In  any  case  the  computational  complexity  and  memory 
storage  requirements  of  this  algorithm  are  excessive.  An 
alternate  was  therefore  examined.  If 


r^  +  r  d  t 


x,  =  r  ,  x,  =  r,  etc. , 


then  a  much  simpler  algorithm  results.  At  a  number  of 


successive  instants  we  have 


(1) 


^2  (1) 


(2)  v-  (2)  +  V,  At 


'll  /  >'A 


and  both  x^  and  x^  may  be  assumed  constant.  If  the  nth 
row  of  the  left  of  (5.28)  is 

v^(n)  V2(n)  +  At  (n-1)  ^ 

T 

then  the  matrix (T^T^^ )  may  be  recursively  formed  as 


(5.27) 


(5.28) 


(T^  T  )  *  (T  /  T  )  +  V  v 

ri  n  n*l  n**x  n  n  • 


(5.29) 
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Although  sufficiently  simple  to  be  feasible  for 
an  air-air  missile,  this  approach  to  range  estimation 
also  failed  for  the  same  reasons: 

(a)  Ill-conditioned  matrix  for  some 
geometries , 

(b)  High  noise-sensitivity  for  short 
averaging  times, 

(c)  Large  errors  in  range-rate  estimates 
for  longer  averaging  times  due  to 
the  Invalid  assumption  that  a™  =0, 

When  failure-mode  (a)  did  not  occur,  and  am  =  0  was 

■^r 

valid,  the  algorithm  yields  estimate  errors  of  the  order 
of  1^, 

No  attempt  at  estimating  a,p  ,  target  acceleration 

■^r 

parallel  to  the  LOS,  by  linear  methods  was  successful. 

Any  algorithm  which  neglected  a™  0  when  it  was  slgnlfl- 

■^r 

cant  operated  unsuccessfully;  either  excessive  noise  sensi¬ 
tivity  became  a  problem,  or  the  neglected  acceleration 
caused  excessive  errors  in  the  estimation  of  Range-Rate, 

In  summary,  except  for  encounter  geometries  which 
yield  ill-conditioned  matrices,  the  matrix-inverse  algorithms 
work  well  with  long-filtering  times  which  eliminate  noise 
effects  if  the  assumption  that  a^  «  0  is  valid.  If  that 
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assumption  is  invalid,  the  matrix-inverse  schemes  fail 
due  to  the  bias  in  range-rate  estimate  caused  by  target 
acceleration,  or  due  to  noise. 

5.3  Observability 

The  failure  of  these  various  adaptive  and  matrix- 
inverse  parameter  identification  methods,  whether  continuous 
or  discrete,  is  fundamentally  due  to  the  existence  of  target 
acceleration  in  the  radial  direction,  parallel  to  the  line  of 
flight.  That  acceleration  component  cannot  be  estimated 
by  linear  techniques  because  it  is  "unobservable"  in  the 
mathematical  sense.  This  section  relates  the  concept  of 
"observability"  to  this  application  and  problem. 

The  state  of  a  linear  system  with  constant  coefficents 
in  the  dynamics  and  constant  coefficients  in  the  output 
relation  can  be  e  ;tlmated  from  knowledge  of  the  output  pro¬ 
vided  the  system  is  mathematically  observable.  For  this 
discussion,  we  view  the  quantities  of  range,  range-rate, 
etc.,  as  states  which  are  to  be  determined  and  consider  the 
signals  of  LOS-rate  and  missile  acceleration  as  time-varying 
parameters.  In  state-variable  notation  a  plant  may  be 
described  as 

Plant  Dynamics;  x  =  Ax  +  Bu 

Output :  y  =  Hx  +  V 

(5.30) 

where  v  is  zero-mean  Gaussian  white  noise,  and  x  is  an 
N-component  vector. 
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The  observability  of  the  system  is  specified  by  a 
relationship  between  the  A  matrix,  specifying  the  dynamics, 
and  the  H  matrix,  relating  the  state  to  the  output.  If 
the  order  of  the  system  is  N,  then  the  matrix 


-  [  H 


T  ,  ,T  „T  * 


H' 


’  (A^)  » 


’(A^)  ] 

t 


(5.31) 


must  have  full  rank  for  the  system  to  be  observable.  Let 
us  now  apply  this  to  one  of  the  possible  formulations  of 
the  problem.  For  the  system  definition,  use  the  accelera¬ 
tion-component  equations  from  Section  2 


accelerations 
parallel  to  the 
LOS 

Perpendicular  to 
the  LOS 


r-ro  +a^  =0 

r  r 


r  o  +  2  f  a  -  a^  +  a^  =  0. 

o  o 


(5.32) 


(5.33) 


In  order  to  cast  this  into  the  standard  form  of  (5.30), 
define  state  variables 


(5.3^) 
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l.e.  the  variables  x  are  the  estimates  of  the  unknown 
parameters.  The  last  two  lines  of  (5«3M  express  the  fact 
that  the  two  components  of  target  acceleration  must  be 
estimated  and  are  assumed  to  be  unknown  constants. 


The  state  variable  expression  of  the  relationships 
among  the  several  variables  may  be  found  from  (5*32)  and 
(5.3^)  as 


(5.35) 

The  first  row  of  (5.35)  expresses  the  differential 

equation  relating  range  and  range-rate  estimates.  The 

second  row  consists  of  (5-3^)  while  the  third  and  fourth 

rows  express  the  assumption  that  the  target's  unknown 

accelerations  a™  and  arp  are  constant.  From  (5.35),  the 
■^r  0 

matrix  A  is 


1 

0 

0 


0  0 


0  o' 
1  0 
0  0 
0  0/ 


and  N=^. 


(5.36) 
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Equation  (5.33)  yields  the  output  relationship 


(5.37) 

comparison  of  the  estimate,  aj^j  ,  with  the  measured 

o 

acceleration  a^^  yields  an  error  whose  minimization  will 
presumably  yield  optimal  estimates  of  the  range,  x^,  and 
range-rate,  X2-  From  (5.37), 


Combining  (5.36)  with  (5.38)  in  the  format  of  (5.35) 
yields  the  observability  matrix 


(5.39) 

,  2 

The  fourth  column  of (5. 39)  is  exactly  equal  to  a 
times  the  second  column,  and  the  rank  is  less  than  the 
order.  This  system  is  unobservable;  experiment  confirmed 
this  result. 
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In  realistic  combat  trajectories  it  is  generally 

erroneous  to  assume  that  a^  ,  target  acceleration  in  the 

r 

range-direction,  is  zero.  If,  however,  that  assumption  is 
valid  for  some  situation,  then  we  can  simplify  (5.35)  by 
omitting  the  third  row  and  column  and  similarly  simplify 
H,  so  that  if  a™  E0,a  priori,  then 


,  =  (-0  -2o  1  ) 


and  0- 


.  *  .21 

-o  -2o  -a  a  I 


-2a  -a  -2o 


In  this  case  rank  equals  order  and  the  system  is 
observable  if  o  ^  0  and  o  0.  This  result  was  experi¬ 
mentally  confirmed. 

Returning  to  (5.35),  assume  that  the  target  accelera¬ 
tion  perpendicular  to  the  LOS,  ,  is  identically  zero, 

A  O 

i.e.,  =  a,p  =  0,  a  priori,  but  the  other  component, 

x^  =  am  ,  is  a  non- zero  constant.  Then,  making  the 
appropriate  changes  in  A  and  K,  we  have 
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and 


0 


32 


-a  -2a  -a 

-2o  -o  -2o* 

0  -2o  0 


.  2 

In  this  case  the  third  column  equals  o  times  the 
first;  therefore  the  system  is  unobservable. 


(5.^1) 


The  result  may  now  be  stated:  If  the  target  has 
significant  acceleration  parallel  to  the  line  of  sight, 
(Sm  y  0),  it  is  not  possible  to  estimate  range  and 

range-rate;  the  system  is  unobservable.  If  the  target 
acceleration  a,-  is  known  to  be  zero,  a  priori;  then  the 
system  is  observable,  and  it  may  be  reduced  to  a  problem 

^  ^  A 

in  the  three  states, x,  =  r,  x-  =  r  and  =  am  .  In  this 
case,  Equation  (5«32)  can  be  exactly  instrumented.  If 

A 

Xn  =  a„  =  constant  is  a  reasonable  assumption,  it  is 
possible  to  reduce  the  problem  from  three  states  to  two 
by  differentiating  the  output  expression  (5.37)  and 

A 

eliminating  =  a^  ,  as  in  (5.16). 

o 

It  is,  however,  an  important  restriction  that  the 
theorems  on  observability  do  not  apply  perfectly  to 
systems  whose  coefficients  are  time-varying.  In  our  case 
the  system  dynamics  are  constant,  while  the  observation 
matrix  is  not.  Intuitively,  it  seemed  possible  that  if 
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the  coefficients  of  the  output  equation  vary  sufficient¬ 
ly  in  the  observation  Interval,  then  a  result  might  be 
obtained  which  is  of  value.  Many  variations  of  the  above 
development,  using  different  models  of  the  system  dynamics 
and  output  equations  were  therefore  examined  to  pursue  this 
question.  The  results  were  disappointing  whenever  a^  ^  0 , 
or  o  was  small.  The  system  Is  unobservable,  despite  the 
uncertainty  of  the  theory  on  this  question,  whenever  a,p  /  0. 
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^.^4  Other  Linear  Systems 

It  is  advantageous  to  unify  the  results  of  the 
above  discussions  on  kinematic  linear  parameter  identifica¬ 
tion  systems  in  order  to  form  a  basis  for  considering  other 
possible  linear  systems.  We  therefore  re-present  the  problem 
from  a  geometrical  viewpoint  and  then  consider  possible  ex¬ 
tensions  beyond  the  scope  of  this  contract. 

Consider  the  equation  of  acceleration  components  per¬ 
pendicular  to  the  LOS,  from  the  kinematics, 


ro  +  2rG  -a^  +aj^_  =  0. 

0  o 


(5.^1) 


Assume  that  the  range  acceleration,  ji,  1b  zero,  and  may 

therefore  be  neglected.  Assume  further  that  a,j,  is  quasi 

o 

constant,  i.e.,  cannot  change  significantly  in  less  than 

1/2  second.  Then  differentiation  of  (5.^1)  with  a  high 

'^1 

pass  filter  such  as  ,  where  is  of  the  order  of 

1/2  second  or  less,  reduces  (5-^l  )  to 


rv,  +  rVo  +  v  =0 

12  o 


(5.^2) 


\n. 

1  "  \t,s  +  1  / 


jT.so  To  \ 

..  X  1  X  I 


T^rnri  ® ;  "s"  ’'o'  tts+t)  ”0  • 

A-  ^ 

If  in  (5.^2)  we  replace  r  and  r  by  the  estimates  r  and  r, 
(5.^2)  may  be  expressed  as  rv^  +  pv^  +  =  e^;  then  the 


(5.^3) 


solution  for  r  and  f  lies  on  a  locus  such  that  e^=0.  In  this 
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i 

I 

I 


i 


I 


case  (5.43)  forms  a  straight  line,  as  shown  In  the 
sketch  below  for^some  particular  instant. 


Fig.  5.1 

Geometrical  Vj.ew  of  Equation  Error. 

The  entire  problem  of  estimation  is  based  on  the  ne- 

cessity  of  finding  another  line  in  the  (r,  r)  space  : 

>  s 

(a)  The  multiple-error  adaptive  process,  by 
using  a  variety  of  values  of  T^  generates  a 
number  of  other  lines  which  are,  more  or  less, 
linearly  independent.  This  approach  may  be 
categorized  as  achieving  distribution  in  the 
T^-la.g  space  at  one  instant  of  true  time. 

(b)  The  various  single-error  adaptive  nets  and  the 
Kalman  net  take  ad\antage  of  the  fact  that,  at 
some  other  instant  than  that  shown  in  Fig.  5.1, 
the  line  of  zero  error  in  the  (r,  r)  space  Is 
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differently  oriented,  if  a  ^  0,and  is  therefore 
linearly  Independent.  These  approaches  thus 
use  the  idea  of  distribution  in  true  time,  t, 
rather  than  in  the  space. 

(c)  The  matrix  inverse  schemes  gather  all  lines 
over  the  computing  Interval,  thus  achieving 
distribution  in  time,  and  then  find  the  solution. 

A  physically  different  approach  to  finding  a  second 

A 

linearly  Independent  line  in  the  (r»r  )  space  is  obviously 
preferable  if  possible.  Stadimetric  ranging  provides  such 
an  approach.  Seekers  vjhich  enable  measurement  of  some 
function  of  the  target  size,  or  area,  such  as  mosaic  IR, 

TV,  area,  or  correlation,  have  this  capability.  In  a 
single  plane  we  have  the  geometrical  relationship 

r8  =  D 

( 

where  r  Range 

6  Angle  Subtended  by  Target,  assumed  a  small  angle 

D  Target  Size 

as  shown  in  the  sketch  below. 
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Assume  that  D  Is  constant;  this  assumption  Is  fair  If 
D  is  the  diameter  of  the  least  circle  which  contains  the 
target  as  viewed  from  the  seeker  in  both  the  maneuver  and 
cross  maneuver  planes.  From  by  differentiating, 

r  0  +  r  e  «  0 

and  y 

re  +  re  *  67. 

^  (5.^5 

''  «■ 

The  choices  of  r  and  r  for  which  e^  =  0  lie  on  a 
straight  line  which  passes  through  the  origin  of  the  (r  ,r ) 
space  at  every  Instant.  Combining  this  with  the  kinematic 


ranging  sketch.  Fig.  5.1,  yields  the  sketch  below. 


Fig.  5.2 

Combining  Kinematic  and  Stadlmetrlc  Ranging. 
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If  we  solve  (5.^^)  and  (5.^5  simultaneously,*  we  get 

A  A 

one  equation  for  r  and  one  for  f.  These  are  independent  and 
their  Joint  solution  may  therefore  be  determined.  It  is 
usually  better  to  solve  the  two  equations  in  a  decoupled  format 
as  the  solution  can  always  be  accomplished  more  rapidly,  more 
simply  and  more  accurately.  A  Least  Magnitude  algorithm  may 
be  used.  They  are  analytically  valid  in  such  cases  and  are 
attractive  for  they  estimate  a  moving  parameter  with  smaller 
lag  than  quadratic  algorithms,  and  thereby  decrease  the  effect 
of  target  accelerations. 

However,-  this  argument  rests  on  the  rather  uncertain  base 
of  the  assumptions.  One  assumption  is  that  the  target  cross¬ 
track  acceleration  (a.^,  )  is  quasi-constant,  i.e.,  changes 

0 

negligibly  within  two  time- constants  of  the  high-pass  filter 
T  s 

i _  ,  In  (5.^2).  The  other  assumption  is  that  r  L  0,  or 

T^s+1 

at  least  the  unknown  portion  (a-,  )  thereof.  In  fact,  in  any 

^r 

maneuver  of  the  target,  and  am  change  continuously  and 

can  change  rapidly.  This  scheme  is  therefore  less  powerful  than 
appears  at  first. 

However,  stadlmetric  data  can  be  valuable  in  reducing 
the  difficulty  of  range  estimation.  Equation  (5.^1)  states 
one  relationship  between  three  of  the  unknowns  in  the  kinematic 
ranging  problem.  Stadlmetric  data  provides  another  relationship, 
so  that  the  effective  total  number  of  unknowns  may  be  reduced 
by  one.  The  possibility  therefore  exists  that  linear  esti¬ 
mation  of  range  and  range  rate  could  be  accomplished  with 

acceptably  small  error,  due  to  the  unobservablllty  of  a,p  ,  by 

r 

combining  kinematic  and  stadlmetric  methods. 
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It  should  be  observed  that  stadlmetrlc  ranging  has  an 
advantage  over  kinematic  ranging.  In  stadlmetry,  the  sub¬ 
tended  angle  6  Is  positive-definite,  and  its  rate  6  Is  zero 
If  and  only  if  r  *=  0;  this  almost  never  occurs.  On  the  other 
hand.  In  kinematic  ranging,  the  objective  of  a  proportional 
navigation  homing  system  Is  to  keep  o  -  0  at  all  times;  conse¬ 
quently  an  ill-conditioned  covariance  matrix  must  result  at 
least  occasionally. 

This  observation  leads  to  describing  the  possible  weak¬ 
ness  in  the  Idea  of  combining  kinematic  and  stadlmetrlc 

At 

ranging  principles.  In  both  schemes,  express  f  as  a  function 
of  r,  and  assume  r  1  0.  Then,  using  (5.^2)  and  (5.^5) 

^  \  "'l  " 

Kinematic:  ^  ^  ^ 

(5.^6) 

Stadlmetrlc:  r  =-(9/9^  r, 

(5.^7) 

under  the  simplifying  assumption  that  =  0. 

o 

As  previously  observed  for  the  stadlmetrlc  method,  in  (5.47) 

9  Is  always  positive,  and  9  is  almost  always  positive. 

•  < 

Neglecting  the  rare  and  anomalous  case  of  9  -  0,  Implying 
increasing  range,  6/0  Is  positive  and  the  solution  for  r 
as  a  function  of  r  passes  through  the  origin  with  a  time- 
varying  negative  slope  which  is  a  quotient  of  non-zero 
finite  quantities.  In  the  kinematic  approach  (5.46),  the 
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coefficient  of  r  will  normally  have  a  negative  value  but  is 
a  quotient  in  which  the  numerator  and  denominator  rapidly 
approach  small  values  and  ultimately  zero,  in  favorable 
geometries,  and  start  and  remain  at  small  values  In  unfavorable 
geometries.  The  slope  of  (5*^6)  and  the  Joint  solution  of 
(5.^6)  and  {^.‘^7),  therefore  are  increasingly  sensitive  to 
instrument  and  computation  noise  as  the  maneuver  becomes 
small.  In  geometric  terms,  the  two  lines  in  the  (r,  r) 
space  of  Figure  5.2  become  parallel  as  the  maneuver  becomes 
small. 

In  consequence,  it  is  evident  that  the  stadlmetrlc/llnear- 


klnem.atic  method  outlined  above  could  be  improved  by  use  of  the 
nonlinear  kinematic  method  described  in  Section 
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6.  TECHNOLOGY  ASSESSMENT 

It  appears  that  linear  methods  of  kinematic  ranging 
can  work  effectively  only  for  favorable  combat  geometries 
and  if  the  sensor  noise  levels  are  quite  low.  This  implies 
high  accuracy  sensors,  and  frequent  ineffectlvlty . 

The  nonlinear  filter  described  in  Section  ^  estimates 
the  target  velocity-vector  components  and  therefore  can  be 
used  to  deduce  the  target  velocity  and  aspect  angle.  The 
filter  also  estimates  the  target's  acceleration  components, 
the  turn-rate, and  its  derivative.  It  therefore  estimates 
the  entire  specification  of  the  target  dynamics,  and  can  be 
used  to  estimate  the  target  future  trajectory  5  and  time  to  go. 
This  may  be  useful  not  only  to  the  air-air  missile  terminal- 
control  guidance  problem  but  also  to  a  variety  of  terminal 
fire  control  systems  in  aircraft  or  helicopters  Instead  of 
missiles.  This  nonlinear  filter  can  also  be  used  to  eliminate 
Impossible  target  accelerations  from  data  provided  by  other 
tracking  systems  such  as  radar  or  laser  rangers. 


Many  of  the  newer  seekers  of  TV  or  IR  types  provide 
some  indication  of  the  target  size  or  shape.  Target  size 
data  can  be  used  as  a  partial  basis  for  a  ranging  system,  as 
outlined  in  Section  5«  In  addition.  Information  on  the 
target  shape  or  its  aspect-angle  can  be  used  to  provide 
clues  enabling  nonlinear  filtering.  For  example,  if  the 
target  is  viewed  from  a  head-on  or  tail-chase  aspect  it  is  clearly 
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difficult  fcr  the  target  to  accelerate  parallel  to  the 
Line  of  Sight;  it  can  accelerate  across  the  line  of  sight 
with  relative  ease. 

Linear  stadimetrlc  ranging  cannot  hy  itself  estimate 
range  and  range-rate,  but  only  their  ratio.  However,  com¬ 
bing  linear  stadlmetry  and  kinematic  schemes  has  the  effect 
of  reducing  the  number  of  unknovm  variables  to  be  determined 
simultaneously.  The  parameter  identification  problem  may  in 
fact  be  reduced  to  a  problem  in  2-space,  without  using  the 
poor  assumptions  that  r,  a,p  ,  and  a„  are  quasi-constant. 

In  this  situation,  two  or  (preferably)  more  independent 
measurement  instants  enable  an  unique  solution.  The  problem 
of  updating  the  hyperplane  of  Pigs.  5.1  or  5*2  to  another 
Instant  is  then  reduced.  It  is  not  eliminated,  therefore 
the  solution  is  not  rigorously  observable;  however,  the 
combination  may  be  useful. 

In  a  review  of  the  present  technology  and  an  estimate  of 
the  trends  and  potentials, the  interaction  of  cost  with  systems 
technology  is  perhaps  most  important  in  the  air-air  missile 
problem.  These  points  stund  out: 

a)  Seekers  of  all  kinds  are  very  expensive;  a 
small  radar  seeker  is  very  costly  in  weight, 
space  and  dollars.  A  TV  seeker  may  be  less 
expensive  in  each  sense,  but  is  not  cheap.  The 
various  sophisticated  IR  seekers  are  cheaper  but 
are  still  quite  costly  compared  to  the  simplest 
IR  seekers. 

(b)  Sensors  such  as  rate  gyros  are  an  order 

of  magnitude  less  costly  than  simple  seekers. 
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c)  The  introduction  of  medium  or  large  scale 
Integrated  digital  and  analog  chips  has  at 
once  improved  the  reliability  of  the  compu¬ 
tation  process  and  reduced  the  cost  by 
orders  of  magnitude. 

Taken  together,  these  elements  suggest  that  cost/ef fectivity 
considerations  lead  to  improving  missile  system  performance 
hy  combining  a  wide  range  of  the  inexpensive  sensors  with  a 
relatively  sophisticated  computation  capability.  This  may 
enable  use  of  a  relatively  unsophisticated  seeker  of  rela¬ 
tively  low  cost. 

The  nonlinear  range/range  rate  estimator  shown  in 
Fig.  4.1  requires  a  considerable  computation  capability,  and 
requires  input  signals  from  a  pair  of  accelerometers.  It 
also  requires  the  missile  airspeed  components  which  may  be 
obtained  from  relatively  inexpensive  sensors  or  may  be 
estimated  from  accelerometer  and  rate-gyro  data.  But  the 
least  sophisticated  seeker  is  sufficient.  Our  approach  was 
to  maximize  the  information  which  could  be  extracted  from 
an  inexpensive  seeker  by  combining  several  Inexpensive  sen¬ 
sors  with  a  relatively  sophisticated  yet  inexpensive 
computer  which  will  soon  be  available. 


86 


Ball  Aaroapace  Company 


REFERENCES 


1.  Eyckhoff,  P.,  "Aspects  of  Process  Parameter  Estimation," 
IEEE  Transactions  on  Automatic  Control,  October,  1963. 

2.  Narendra,  T.  S.,  and  L.  E.  McBride, "Multi-Parameter 
Self-Optimizing  System  Using  Correlation  Techniques," 
IEEE  Transactions  on  Automatic  Control,  January,  I962. 

3.  Young,  P.  C.,  "The  Use  of  Linear  Regression  and  Related 
Procedures  for  the  Identification  of  Dynamic  Processes," 
IEEE  Proceedings  of  the  Seventh  Symposium  on  Adaptive 
Processes,  UCLA,  December,  1968. 

Lion,  P.  M. ,  "Rapid  Identification  of  Linear  and 
Nonlinear  Systems,"  AIAA  Journal,  Vol.  5,  No.  10, 
October,  I967. 

5.  Bell  Aerospace  Company  Report  9500-920183,  "Adaptive 
Modeling  and  Control,"  F.  D.  Powell,  June,  1970, 

Final  Report  on  Office  of  Naval  Research 

Contract  Nonr-486^ ( 00 ) ,  NR0^9-207. 

6.  NASA  Report  68-11-7,  "Use  of  Integral  Transforms  in 
Parameter  Estimation."  H.  C.  Lessing  and  D.  F.  Crane, 

1968. 

7.  Morrison,  N.,  "Introduction  to  Sequential  Smoothing 
and  Prediction,"  McGraw-Hill  Book  Company,  New  York 

1969. 

8.  Kalman,  R.  E. ,  "A  New  Approach  to  Linear  Filtering 

and  Prediction  Problems,  "  Transactions  of  the  A.S.M.E., 
Journal  of  Basic  Engineering,  Vol.  82,  PP  35-^5, 

March,  I960. 


Bell  Aeroepaee  Company 


88 


B«ll  A«ro«paoft  Company 


We  also  demonstrated  that  the  acceleration  compoiient 
equations  parallel  and  perpendicular  to  the  LOS  are, 
respectively : 


.  a 

r-ra  -a^  +  aj^,  "0 


and 


rd  2ra  -  a^  +  -  0 

o  o 


The  target  and  missile  velocity  and  acceleration 
components  parallel  and  perpendicular  to  the  LOS  are: 

Velocities : 


Parallel;  cos 

r 


(A. 3) 


(A. 4) 


Perpendicular;  V„  »  V„  sin 


Accelerations ; 


Parallel;  a^  -  cos  -  V^yt 
r 

r 

Perpendicular:  a,j,  -  sin  5ij.  +  V^  y^  cos 

0 


a 


M 


sin  +  Vj^Ym  cos 


(A. 5) 


(A. 6) 


where 
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Angles : 


(A. 7) 


The  principle  of  this  approach  is  simply  to  solve 
A, 3  for  r  as 

.  2 

r  *  ro  +  aj^  -  a,p  (A.  8) 

r  r  • 

Integration  of  r  will  yield  r  and  r  if  the  terms  on  the 

right  of  (A. 8)  are  available. 

.  2 

The  term  ro  can  be  computed  easily  while  a.^  ,  the 
missile  acceleration  parallel  to  the  LOS,  can  be  measured 


directly  The  problem  therefore  is  nov;  to  estirr.ate 


in  some  way.  It  is  assumed  that  initial  values  of  r  and  r 
are  provided. 

It  was  shov/n  that  the  target  acceleration  in  the  range 


direction,  a„  ,  can  be  determined  from  available  data.  The 
^r 

procedure  is  first  to  calculate  V^,  the  target  airspeed,  by  use 
of  (A. 1)  and  (A. 2)  as 


(f  +  cos  ^j^)^+(ro  +  Vj^,  sin  = 

cos  sin  =  V^.  ^ 

It  was  stated  in  Section  (^.1)  that  this  expression 
may  be  differentiated  to  yield  a^^.  +  Y„  a„  =  V_  V„; 

This  is  demonstrated  below.  Therefore,  let  us  differentiate 
(A. 9)  with  respect  to  time,  yielding 

2{(f+Vj^^  cos  5j^)Cr+(V^^  cos  sin  sin 

+  (ro+Vp^  sin  Sj^^)[rd+ro  +(Vj^  sin  cos  C|y,)  - 

=  2  V^V^. 
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Substitute  from  (A. 5)  the  simplifying  notation  for  the 


acceleration  components  of  the  missile,  a,,  and  a 


parallel  and  perpendicular  to  the  LOS,  respectively;  these 
accelerations  are  gathered  in  parenthesis  within  the  same 
brackets  in  the  expression  above  [(  )]. 


This  yields 


(r  +  Vj,  cos  Cf^j)(r  +  a^^  +  sin 

r 

+  (ro  +  Vj^,  sin  Cj;|)(rd  +  ro  +  aj^^  -  Vj^o  cos  = 


(A. 1C) 


Now  solve  (A.  2)  for  V,_,  sin  ,  multiply  by  o  and  replace 
the  term  (Vj_.o  sin  j  on  the  first  line  of  (A. 10)  by  th.e 
term  thus  formed.  Similarly,  solve  (A.l)  for  V,,  cos 
multiply  by  a  and  eliminate  the  term  (v,.  o  cos  ^„)  from  the 

j*l  *  * 

second  line  of  (A. 10).  These  substitutions  yield 

.  ...  2 

(r  +  V,.  CCS  Cf,j)(r  -  ro  +  +aVrp  sin  £;,p,)  + 

‘‘r 

(ra  +  V,,  sin  5„)(rd  +  2ro  +  a..  -oVr..  cos  f  )  = 

I'i  1'!  i  1 

V  V 
''  T  T 

Keeping  the  right  side  of  this  expression,  gather  the 

coefficients  of  on  the  left;  this  reduces  to 

.  .  2 

(r+V,,  cos  C.J(r  -ro  +  a.,  )  +  (ro  +  V.,  sin  C,.) 

I'I  ^i  I'l^  I'I  I'l 

(rd  +2ra  +  a,,  ) +o  {(r+v..  cos  C,.)(Vrr-  sin  Crr)  - 

i'l  I'l  I'l  I  i 

o 

-(ro  +  sin  cos  C^))  =  V^V^. 


(A. 11) 
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The  term  in  curly  braces  {  }  in  (A. 11)  is  identically 
zero.  This  follows  from  (A.l),  r  +  cos  ^  cos  C,p, 

while  from  (A.  2),  ro  +  Vj^  sin  sin  so  that  the 

brace  has  the  value  {V^  cos  sin  sin  cos 

As  a  consequence,  (A. 11)  takes  the  simpler  form 

(f+V^  cos  Cf^)(r-ro  +  ajj^  )  +  (ro+Vj^^  sin  Cj^,) 

r 

(ro  +  2ra  +  a„  )  =  V„V  . 
r 

By  use  of  (A.l)  through  (A. 6),  this  may  be  expressed 

in  the  intuitively  satisfying  form  V  a^  +  V  a^  =  V^V^pj 

■^r  ■^r  "^o  "^o 

which  is  the  form  used  in  the  discussion  in  Section 
where  we  use  the  simplifying  riotatlon,  from  (A.  5), 

cos 

r 

Vt  =  sin 

rr 

to  express  the  target  velocity  componen'  3  parallel  and 
perpendicular  to  the  LOS. 


The  exact  solution  for  a,j,  is  thus 

■^r 

•  V  a 
■  V„ 


and 


•  V  a 

Vt  -  V  X 

r  =  I’D  -  aj,^^  + 


while  a„ 


— ^  =  Y  if  V  “  0 
Y  ’T  ''t 


0. 

(A. 12) 

(A. 13) 


(A.  in) 

(A. 15) 
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Usually  Is  small  and  can  be  neglected,  and  all  other 
terms  on  the  right  of  (A.  15)  can  be  estlmat^ed. 

The  mechanization  of  this  algorithm  is  discussed  in 
Section  (A. 2),  following. 


A.  2  Mechanization  of  the  Slg,..ais  for  the 
_ Ideal  Algorithm _ 

The  range  acceleration  is  given  by  (A. 15);  it  may  be 

.  .  2 

integrated  to  yield  r  and  r.  The  term  ro  can  then  be 

formed,  while  the  missile  acceleration  component,  aj^^j  ,  can 

r 

be  directly  sensed. 


The  mechanization  of  (A. 15)  requires  that  missile 

airspeed  Vj^,  the  missile  acceleration  components  aj^  and 

i  aj^  ,  and  the  angle  be  available  as  signals  or  be  estl- 

o 

I  mated.  It  will  now  be  shown  that  this  can  be  accomplished. 


(a)  Acceleration  Components 

The  accelerometer  components,  a^^  ,  and  a^^  ,  are  a 

r  o 

special  case;  two  easy  solutions  are  available; 

(1)  Mount  the  accelerometers  on  the  seeker; 

then  one  yields  and  the  other  yields 

r  o 

directly . 

(2)  Mount  the  accelerometers  on  the  missile 

body;  then  their  signals  are  aj^^  ,  and  aj,,  ,  fore  - 

X  y 

aft  and  cross-body,  respectively;  and  resolution 


of  the  signals  from  the  accelerometer  through 
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the  "look-angle,”  X,  yields 


/  \ 

[  cos  X  -sin  X  1 

ixl 

^sln  X  cos  X  j 

As  the  "look-angle,”  X,  between  the  missile 
centerline  and  the  LOS  can  easily  be  made 
physically  available  in  a  seeker  gimbal,  the 
resolution  of  the  acceleration  components  is 
easily  accomplished. 

(b )  Velocity  Estimation 

If  the  missile  body  pitch  rate,  q,  or  equivalent, 
is  an  available  signal,  and  further,  if  the  cross¬ 
body  acceleration  of  the  missile,  a  .  ,  is  also 

'■'‘y 

available  as  a  signal,  it  is  relatively  easy  to 
estimate  missile  velocity,  . 

The  pertinent  equations  of  motion  of  the  missile 
are , 


a|^ij^  =  Nci  +  G6 

where  a  angle  of  attack 

aj^^  cross-body  acceleration 

‘y 

6  control  displacement 

q  body  pitch  rate 

airspeed 


9k 


A. 16) 


A.  17) 
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N,  G  cross-body  accelerations  per 

radian,  due  to  angJ.e  of  attack 
and  control  displacements, 
respectively 

Assuming  quasi-constant  conditions,  the  transfer 
function  defining  acceleration  is 

_  N  a  +  Gs6 
-  s  +  N/V^. 

The  acceleration  due  to  control  displacement  may 
usually  be  neglected  so  that 

'  N/V, 

=  V,,  1 

■y  'M 

s  V  _ 

K.  Ts+1 


M. 


s  +  N/V 


M  i  s+N/Vj^^l 


(A. 16)  I 


(A. 19) 


where  T  =  •  Then  Vj,^  may  be  estimated  by  the  algorithm: 


or 


Vf,  =  k  a,.,  (a,.  -  V,,  q) 


V,,  =  k  (  I  a..  I  -  V..  I  a  I  ) 

I'l  •  1<1  '  Iv  '  •  ' 


1>. 


where  k  >  0  has  the  appropriate  value. 

Somewhat  superior  algorithms  which  will  lag  less 
are  achieved  by  biasslng  the  algorithrus  given  by  tlie  missile 
longitudinal  acceleration  component  a,,  . 


(A. 20)  i 

o 

(A. 21)  I 
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Algorithm  (A. 21)  Is  obviously  simple  to  instrument, 
requiring  only  one  multiplication. 

(c)  Estimation  of  the  Angle 

The  angle  5jv]  is  defined  as 
But  as  Yjvi  is  the  centerline 

attitude,  then  =  e^^-a-a,  and  as  the 
angle  from  missile  centerline  to  the  seeker  LOS  is 

X  .  <7  -  0,,, 
we  have  =  -^  “  a 

where  a  is  the  angle  of  attack  between  the  missile 
centerline  and  the  velocity  vector.  But  this  can 
easily  be  estimated  as 


(A. 22) 


Ts  +  1 


vjhere  T  was  defined  above, 


(A. 23) 


The  missile  velocity  components  Vj^  and  Vj^ 

r  ‘  o 

can  then  be  computed  as  *  V.  cos  and 
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APPENDIX  B 

Error  Analysis  Coefficients 

The  numerical  data  presented  in  this  appendix  represent 
the  time  varying  cpefflclents  of  the  error  equation.  As  dis¬ 
cussed  in  Chapter  ,  the  error  equation  Is  of  the  form 

X  +  a, X  +  a-x  *  )  b .N. , 

1  2  1  1* 

and  the  data  here  represent  the  coefficients  and  tabulated 
for  every  0.25  seconds.  There  are  six  tables,  representing  a 
selection  of  six  choices  of  initial  geometry  configuration.  They 
are  ordered,  and  correspond  respectively  to  Figures  il.3  through 
^.8  incluslvej  however,  each  graph  is  labelled  to  show  the  initial 
geometry.  These  data  permit  more  general  stochastic  analyses  of 
the  effects  on  the  estimates  due  to  errors  In  the  ,  'eral  signals 
which  are  used  in  the  algorithm.  The  coefficients  ^  and  are: 

a^^  =  2o  tan  tan  e,p  ,  and 

•  ^  *  •  .. 

32  =  -0  +  TrpO  +  0  tan 

The  coefficients  b^  and  error  sources  are: 

b.  *  -1  N,  a„  ,  accelllLOS  error 

1  1  Mj,* 

b2  *  sec  Crp  N2  V,p,  target  airspeed  change 

A. 

b^  “  -Y.J,  sec  sin  (Yf,^-y,p)  Vj^,  missile  airspeed  error 
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b|j  »  -rVip  -2r  tan  Cfj.  +  2ra 

bg  -  -r  tan  5.p 
by  =  -tan  Cy, 


N^j  o,  LOS  rate  error 

5jv|»  error  of  estimating  a 
Ng  o  LOS  acceleration  error 
Ny  aj^  ,  accel.  -L  los,  error 


All  angles  are  In  radian  measure  unless  specifically  stated 
otherwise . 

The  data  of  this  appendix  are  presented  in  a  conventional 
"exponential  format"  for  digital  output  data.  Thus  the  tabulated 
item  0.^22E  01  Is  interpreted  as  (0.^22)10"’’^  *  ^.22,  while  -0.2it^E-01 
is  interpreted  as  (-0.24^1)10”^  »-0.024U. 
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The  matrix  is  thus  a  "continuously  updated  observa¬ 
tion  matrix"  and  the  vector  is  a  "series  of  output 

observations".  The  expression  (S.21)  is  valid  provided 
the  noise  v  is  Independent  from  computing  instant  to 
instant  and  has  zero  mean.  The  construction  of  the 
matrix  T^  is  often  a  formidable  task  as  it  requires 
keeping  in  storage  a  large  amount  of  past  data  K(t^)  so 
that  it  can  be  updated  by  the  state  transition  matrix 
<I>(t.jtj)  which  is  at  times  also  difficult  tc  obtain. 

1  -j 

Apply  this  theory  to,  typically,  the  rcd\iced 
dimensionality  problem  defined  by  (5.16)  and  (5-17); 
aSi;uiue  tliat  r  iis  eonbtanl,  t'nu.s  X-  is  constant  afbei’ 
convergence.  But  =  r  and  is  not  constant.  Now  we 
have,  at  some  specific  instant,  and  with  only  one  value 
of  1  =  1 , 


v^(t^)x^  +  v^Ct^)/.^  +  y^  = 

At  the  next  instant  we  have 


(5.2^) 


T  (  \  ^  V 

1  "'2'"1 


^'2~  ^5' 


But  while  the  value  of  is  Ideally  the  same  at 
these  two  in.stants,  the  value  of  chanp'es  continuously, 
for  -  r,  and  r  is  the  integral  of  r.  V/e  may  modify 
(5 .2^)50  that  is  the  same  for  the  two  Instants.  If  At 
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i' s  the  time  between  samples,  vie  may  rewrite  (5.24)  as 
v^(t^)  [xAx..^')  +  x„  At)  +  v„  {t^)  =■•  Cp  , 

which  may  be  rearraiiged  as 

Vi(Xi)  x^(t2)  +  (v^  (t^)  +  (t^)  At)  “  ^'5 

Now  Eqns .  (5.25)  and  (5-26)  represent  two  ineasure.mer.t.s 
referred  to  the  same  instant,  and  may  be  solved  Jointly, 
liore  generally,  a  larger  sequence  of  such  instants  p-.ay  be 
gathered  in  the  more  general  form  of  (5.23). 

This  method  cxhioited  several  modes  of  failure. 
Whenever  any  two  of; 


the  missile  velocity  vector 
the  target  velocity  vector 


the  LOT 


become  collnear,  or  c  -•  0 ,  the  n.atrix  (T  'T  )  became  ill- 

n  n 

conditioned  or-  indeterminate,  yielding  inaccurate  and  noise- 
sensitive  solutions.  In  geometries  when  this  mode  of  failure 
does  not  occur,  this  algcrlthn!  has  a  different  n:cde  of 
failure.  If  the  durario:-|  of  t.he  st  rn*?iT'or'y 

1/2  second,  the  neglected  target  acceleration  parallel  to 
the  LOS,  ,  causes  the  estimate  of  r  to  be  tc  5C5  in 

1  ^  ^  I 

V 

error  within  1  to  2  seconds,  as  the  assumption  that  f  =  0  Is 
too  Incorrect.  On  the  ath;-r  hand,  a  shorter  duration  memory 
results  in  high  noise  sensitivity.  but  triese  algoiMtliP.s  v.'ork 
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well  (Range-Rate  errors  of  1%)  if  a^p  =  0. 

r 

In  any  case  the  computational  complexity  and  memory 
storage  reo.uiremcnts  of  this  algorithm  are  excessive.  An 
alternate  was  therefore  examined.  If 


d  t 


and  X,  =  r  ,  X.  =  r,  etc., 

1  o  ’  ^  *  * 

then  a  much  simpler  algorithm  results.  At  a  number  of 
successive  instants  we  have 


and  both  x^  and  x^  may  be  assumed  constant.  If  the  nth 
row  of  th.e  left  of  (5-28)  is 


(5.28) 
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Although  sufficiently  simple  to  be  feasible  for 
an  air-air  missile,  this  approach  to  range  estimation 
also  failed  for  the  same  reasons: 

(a)  Ill-conditioned  matrix  for  some 
geometries , 

(b)  High  noise-sensitivity  for  short 
averaging  times, 

(c)  Large  errors  in  range-rate  estimates 
for  longer  averaging  times  due  to 
the  invalid  assumption  that  =  0, 

When  failure-mode  (a)  did  not  occur,  and  arp  *  0  vias 

r 

valid,  the  algorithm  yields  estimate  errors  of  the  order 
of  1%. 

No  attempt  at  estimating  a^,  ,  target  acceleration 

r 

parallel  to  the  LOS,  by  linear  methods  was  successful. 

Any  algorithm  which  neglected  /  0  when  it  was  slgnlfi- 

cant  operated  unsuccessfully;  either  excessive  noise  sensi¬ 
tivity  became  a  problem,  or  the  neglected  acceleration 
caused  excessive  errors  in  the  estimation  of  Range-Rate, 

In  summary,  except  for  encounter  geometries  which 
yield  ill-conditioned  matrices,  the  matrix-inverse  algorithm 
work  well  W'lth  long-filtering  times  which  eliminate  noise 
effects  if  the  assumption  that  =  0  is  valid.  If  that 


>-i  •> 
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assumption  is  invalid,  the  matrix-inverse  schemes  fail 
due  to  the  bias  in  range-rate  estimate  caused  by  target 
acceleration,  or  due  tc  noise. 


5.3  Observability 


The  failure  of  these  various  adaptive  and  matrix- 
inverse  parameter  identification  methods,  whether  continuous 
or  discrete,  is  rundamcntally  due  to  the  existence  of  target 
acceleration  in  the  radial  direction,  parallel  to  the  line  of 
flight.  That  acceleration  component  cannot  be  estimated 
by  linear  techniques  because  it  is  "unobservable"  in  the 
mathematical  sense.  This  section  relates  the  concept  of 


"  ub  h  f  r  v  a’o  lilt  y ' 


p  X  j.  V  a  V  JL  Vi  1  aiiu  pxvwxUi 


The  state  of  a  linear  system  v;ith  constant  coefficents 
in  the  dynamics  and  constant  coefficients  in  the  output 
relation  can  be  estimated  from  knowledge  of  the  output  pro¬ 
vided  the  system,  is  mathemiatlcally  observable.  For  this 
discussion,  we  view  the  quantities  of  range,  range-rate, 
etc.,  as  states  which  are  to  be  determiined  and  consider  the 
signals  of  LOS-rate  and  missile  acceleration  as  time-varying 
param.eters.  In  state-variable  notation  a  plant  may  be 
described  as 

Plant  Dynamics;  x  =  Ax  +  bu 

Output ;  y  =  Hx  +  V 

(5. 

where  v  is  zero-mean  Gaussian  white  noise,  and  x  is  an 
IJ-ccrr.ponent  vector. 
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The  observability  of  the  system  is  specified  by  a 
relatlonsliip  between  the  A  matrix,  specifying  the  dynamics, 
and  the  H  matrix,  relating  the  state  to  the  outpnt.  If 
the  order  of  the  system  is  N,  then  the  matrix 

'  '  ,  '  '  n-l 

rp  rnm  m2m 

0,,  =  [  H*  '  •  (A^)  '  .  .  .  '  (A  ;  ] 

»  I  I  t 


(3.3:) 


must  have  full  rank  for  the  system  to  be  cbsei'vable.  Let 
us  now  apply  this  to  one  of  the  pcssitlc  formulations  of 
Che  problemi.  For  the  system  definitior. ,  use  the  accelera¬ 
tion-component  equations  from  Section  2 


accelerations 


parallel  to  the 


r  -  r  o  -  a,,  =0 


Ferpendlcular  to  r  o  +  2  r  c  -  arp  +  a., 


M 

0  0 


the  L0.S 


(5.32) 


(0 


In  oraer  to  ca.sl  this  into  the  standard  fcrmi  of  (5-30), 
define  state  variables 


(5.  3'0 


71 


Bell  Aerospace  Company 


i.e.  the  variables  x  are  the  estimates  of  the  unknown 
paraiTieters .  The  last  two  lines  of  (5.3^)  express  the  fact 
that  the  two  components  of  target  acceleration  must  be 
estimated  and  are  assumed  to  be  unknown  constants. 


The  state  variable  expression  of  the  re latlonslilps 
among  the  several  variables  may  be  found  from  (5.32)  and 
(5.3^)  as 


The  first  row  of  (5.35,1  expi’esses  the  differential 
equation  relating  range  and  range-rate  estimates.  The 
second  row  consists  of  (5-3^)  while  the  th.ird  and  fourth 
rows  express  the  assumption  that  the  target's  unknown 
accelerations  a,,  and  a,p  are  constant.  From  (5.35),  the 
matrix  A  Is 


and  N=^. 


(5.36) 
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acceleration  aj^  yields  an  error  whose  minimization  will 
^  o 


presumably  yield  optimal  estimates  of  the  range,  x^,  and 
range-rate,  X2-  (5-37), 

H  =  (  -0  -2o  0  1  ). 


(5.38) 


Combining  (5-36)  with  (5-38)  In  the  format  of  (5.35) 
yields  the  observability  matrix 


-a 

,  3 

-?o 

1 

Q  • 

-2a'  \ 

•  0 

(M 

-d 

.  S 

-2o 

-0  0 

0 

~2h 

-a 

m 

•  D 

CM 

1 

1 

0 

0 

0  / 

\  / 


(5.39) 

•  ^ 

The  fourth  column  of (5- 39)  is  exactly  equal  to  o 
times  the  second  column,  and  the  rank  Is  less  than  the 
order.  This  system  Is  unobservable;  experiment  confirmed 
this  result. 


4 


i 
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In  realistic  combat  trajectories  it  is  generally 
erroneous  to  assume  that  a„  ,  target  acceleration  in  the 
range-direction,  is  zero.  If,  however,  that  assumption  i 
valid  for  some  situation,  then  we  can  simplify  (5.35)  by 
omitting  the  third  row  and  column  and  similarly  simplify 
H,  so  that  if  a^  =0,a  priori,  then 

o\ 


A  ^  — 

^  JL 


r 

0  1 

.  2 

a 


0 

0 


*  ^31 


(-0  -2o 


1  ) 


and  0 


31 


-0 

-2a 

\  ^ 


.3  .2 

-2o  -c  0 

.  3 

-d  -2o 

0  0 


In  this  case  rank  equals  orutrr  and  the  system  is 
observable  If  d  0  and  o  0.  This  result  was  experi¬ 
mentally  confirmed. 


Returning  to  (5>35)s  assume  that  the  target  accelera 

tion  perpendicular  to  the  LOS,  a^p  ,  Is  identically  zero, 

i.e.,  =  0,  a  priori,  tut  the  other  component, 

=  a^p  ,  is  a  non-zero  constant.  Then,  making  the 
■^r 

appropriate  chajiges  j.ii  k  and  }i,  we  have 


32 


0 

.  2 

0 


1 

0 

0 


1 

c 


‘32 


(  -d  -20  0) 


I 

I 

L 
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and 


-2o 

II 

f\J 

o 

—  C.0 

-d 

1  ° 

• 

-  2g 

In  this 

cas e  the 

third 

first;  there I'cre  the  syste.r. 


-o  a 


0 


/ 

.  2 

column  equals  o  tlmies  the 
Is  unobservat le . 


The  result  may  now  be  stated:  If  ti.e  larfet  has 


significant  acoc  Icrat  ion.  parallel  to  the  line  of  sight, 
(a^  ¥  0),  it  is  not  possible  to  estimate  range  arid 


r 

range-rate;  the  systein  is 

acceleration  ar*.  is  kncwri 
i 

systemi  is  observable,  and 

in  the  three  states, = 

case,  Kquatlor;  (5.32)  can 

X ,,  =  =  constant  is  a 

i 

o 

possible  tc  reduce  the  pr 


unobservable.  If  the  target 
tc  be  zero,  a  priori;  thei:  ths 
it  !;.a,v  be  reducec  to  a  protiemi 

^  ^  /N 

r,  -  1'  and  x ,,  =  a^.  •  Ir  th.i 

c:  4  1^ 

O 

be  exactly  Inst  run, ented  .  If 
reasonable  assumption,  it  is 
cblem  from  tb.ree  states  to  tv.o 


by  different  i.at  ing  the  output  expression  (5.37)  and 

eliiTiinati  f.g  X;j  =  ,  as  In  (5.1b). 

'  a 

It  is,  however,  aii  importarit  restrlcticn  that  tb.e 


theorems  on  observability  do  not  apply  perfectly  to 
systems  whose  coefficients  are  time-varying.  In  our  cas 


the  system  dynamics  are  constant,  while  the  observaticr; 
matrix  is  not.  Intuitively,  it  seemed  possible  tiiat  if 
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the  coefficients  of  the  ovitput  equation  vary  sufficient¬ 
ly  in  the  observation  interval,  then  a  result  might  be 
obtained  v.’hich  is  of  value.  Hany  variations  of  the  above 
development,  using  different  models  of  the  system  dynamics 
and  output  equations  were  therefore  examined  to  pursue  this 
question.  The  results  were  disappointing  whenever  arp  ^  0 , 
or  0  was  small.  The  system  is  unobservable,  despite  the 

uncertainty  cf  the  theory  on  this  question,  whenever  a,,.  /  f . 

'  r 
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5.^  Other  Linear  Systems 

It  is  advantageous  to  unify  the  results  of  the 
above  discussions  ori  kinematic  linear  parameter  identifica¬ 
tion  systems  in  order  to  form  a  basis  for  considering  other 
possible  linear  systems.  V/e  therefore  re-present  the  problem 
from  a  gecmetrjcal  viewpoint  and  then  consider  possible  ex¬ 
tensions  beyond  the  scope  of  this  contract. 

Consider  the  equation  oj‘  accelerallcii  components  per¬ 
pendicular  to  the  LOS,  from  tiie  kinemiatics, 


r  0  +  2  r  G 


Assume  that  the  range  acceleration,  ,  is  zero,  and  may 

therefore  be  neglected.  Assume  furtlier  that  a,.  is  quas 

0 

constant,  i.e.,  cannet  change  sign! f leant ly  in  less  than 

1/2  second.  Then  differentiation  of  )  with  a  high 

'^1 

pass  filter  such  as  »  v.-here  is  of  the  order  of 

1/2  second  or  less,  reduces  (5.^1)  to 


r  V  -  +  >-  V  T  +  V  =  0 

i  ?  c 


where  v 


i  S  . 1 ) 


'  j  *  ^  j 


solution  for  i'  and  r  lies  on  a  locus  such  that  e—0.  In  thi 


_  pT.si  ^  T,3 

\  „  _ ! 

''i^  \a 

1  T" 

V 

T^s+1  /  ’  2 

V^i^+l  (T^s+ 

1)^’  'o  ( 

T^s  +  l]  "o' 

1 

)  we  replace 

r  and  r  by  t'n.e 

estimate 

s  r  and  r , 

8 

be  expressed 

as  r’/^  +  + 

V.  =  e.; 

then  the 

{ 

h  .  A  ' 
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case  (S^3)  forms  a  straight  line,  as  shown  in  the 
sketch  below  for^sorr.e  particular  instant. 


Fig.  5.1 

Geometrical  View  of  Equation  Error. 

The  entire  problem  of  estimation  is  based  on  the  ne- 

cessity  of  finding  another  line  in  the  (r,  r)  space: 

••  >  6 

(a)  The  multiple-error  adaptive  process,  by 
using  a  variety  of  values  of  T^  generates  a 
number  of  other  lines  which  are,  more  or  less, 
linearly  independent.  This  approach  may  be 
categorized  as  achieving  distribution  in  the 
T^-lag  space  at  one  instant  of  true  time. 

(b)  The  various  single-error  adaptive  nets  and  the 
Kalman  net  take  advantage  of  the  fact  that,  at 
some  other  Instant  than  th:at  shov.n  an  Fig.  5.1, 
the  line  of  zero  error  in  the  (r,  r;  space  is 
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differeritly  oriented,  if'  o  ^  C,and  Is  therefore 
linearly  independent.  These  approaches  thus 
use  the  idea  cf  distribution  in  true  tine,  t, 
rather  than  in  the  space. 

(c)  The  matrix  inverse  schemes  gather  all  Ij.nes 
over  the  conputing  interval,  thus  achieving 
distribution  in  time,  and  then  find  the  solution. 

A  physically  different  approach  to  finding  a  second 
linearly  independent  line  in  the  (r^r  )  f;pace  is  obviously 
preferable  if  possible.  Stadlm.etrlc  rar^ging  provides  such 
an  approach.  Seekers  v;}.ioh  enable  mear.urement  of  sc.rr.e 
function  cf  the  target  size,  or  area,  such  as  mcsaic  IR, 
TV,  area,  or  correlation,  have  this  capability.  In  a 
single  plane  v;e  have  the  geometrical  relationship 

r&  =  D 


where  r  Range 

0  Angle  .Sutteridf'd  by  Target,  assumed 

D  Target  .Size 


snail  atifl 


as  shovv'n  in  th.e  sketch  beii..v.' 
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Assur.e  that  D  Is  constant;  this  assumption  is  fair  if 
D  is  the  diameter  of  the  least  circle  which  contains  the 
target  as  viewed  from  the  seeker  in  both  the  maneuver  and 
cross  maneuver  planes.  From  (5.^^),  by  differentiating, 


r  e  +  r  6  =  0 


r  6  +  r  e  =  e- 


(5.^5) 


The  choices  of  r  and  r  for  which  e-  =  0  lie  on  a 

( 

straight  line  which  passes  tlirough  the  origin  of  the  (r  ,r  ) 
space  at  every  Irnstant.  Combining  this  v.’lth  the  kinematic 
ranging  sketch,  fIk.  b.lj  yields  the  sketch  below. 


-hi 


,T  o  "I  ^ 

Solution 


Fig. 

Combining  Kinematic  and  Stadimotric  Ranging. 
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II'  we  solve  (5-^^)  and  (5-^5  simultaneously,-  we  f,et 

/V 

• 

one  equation  for  r  and  one  for  r.  These  are  Independent  and 
their  Joint  solution  may  therefor-e  be  determined.  It  is 
usually  better  to  solve  the  two  equations  in  a  decoupled  format 
as  the  solution  can  alv.-ays  be  accomplished  more  rapidly,  more 
simply  and  more  accurately.  A  Least  Magnitude  algorithm  mjay 
be  used.  They  are  analytically  valid  in  such  cases  and  are 
attractive  for  they  estimate  a  moving  parameter  with  sm;aller 
lag  than  quadratic  algorithm;G ,  and  tliereby  decrease  the  effect 
of  target  accelerations. 

Hov.’ever,  this  argumier.t  rests  on  the  rather  uncertain  base 
of  the  assumptions.  One  assumption  is  that  the  target  ci’oss- 
track  acceleration  (a,p  )  is  quasi-constant,  i.e.,  changes 
negligibly  within  two  time- constants  of  the  high-pass  filter 
T  s 

i  ,  in  (5*^2).  The  other  assumjpticn  is  that  r  2  C ,  or 

T^s-t-1 

at  least  the  unknown  portion  (am  )  thereof.  In  fact,  in  any 

^r 

maneuver  of  the  target,  and  a,p  change  continuously  and 

can  change  rapidly.  This  schem^e  is  therefore  less  pcv/erful  than 
appears  at  first. 

Kev.-ever,  stadimetrlc  data  can  be  valuable  in  reducing 
the  difficulty  of  range  estlm.atlcn.  tquation  (5-^l)  states 
one  re j.allonohip  between  t.hree  of  the  unknowns  in  the  kinematic 
ranging  problem.  Stadimetrlc  data  provides  another  relationship, 
so  that  the  effective  total  nurndver  of  unknowns  miay  be  r-educed 
by  one.  The  possibility  therefore  exists  that  linear  esti¬ 
mation  of  range  and  range  rate  could  be  accomplished  with 
acceptably  small  error,  due  to  the  unobservabillty  c^f  a^j,  ,  by 
comibinlng  klnem:atlc  and  stadimetrlc  methods. 
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It  should  be  observed  that  stadlmetric  ranging  has  an 
advantage  over  kinematic  ranging.  In  stadimetry,  the  sub¬ 
tended  angle  e  is  positive-definite,  and  its  rate  6  is  zero 
if  and  only  if  r  -  0;  this  almost  never  occurs.  On  the  other 
hand,  in  kinematic  ranging,  the  objective  of  a  proportional 
navigation  homing  system  is  to  keep  o  i  0  at  all  times;  conse¬ 
quently  an  ill-conditioned  covariance  matrix  must  result  at 
least  occasionally. 

This  observation  leads  to  describing  the  possible  weak¬ 
ness  in  the  idea  of  combining  kinematic  and  stadlmetric 
rai'jging  principles.  In  both  schemes,  express  r  as  a  funclicn 
of  r,  and  assume  v  ~  0.  Then,  using  (5.^2)  and 


Kinemat 1 c ; 


Stadlmetric  :  r  = -(6/6)  r^ 

( '3 . 

under  the  simplifying  assumiption  that  a.-,  =  0. 

0 

As  previously  observed  for  the  stadimetrlc  method,  Ir;  (3. A 
6  lo  alv;ays  positive,  and  6  Is  almost  always  positive. 
K'eglccting  the  rare  and  anomialous  case  of  6  -  0,  implying 
ir'.creaslng  range,  6/6  i.s  positive  and  the  solution  for  f 
as  a  function  of  r  passes  througli  the  origin  v.j.th  a  timie- 
varylng  negative  slope  which  is  a  quotient  ol'  non-zero 
finite  quantities.  In  the  kinematic  approach  (3.A6),  the 
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coefficient  of  r  v.ill  normally  have  a  negative  value  but  is 
a  quotient  in  which  the  numerator  and  denominator  rapidly 
approach  small  values  and  ultimately  zero,  in  favorable 
geometries,  and  start  and  i-emain  at  small  values  in  unfavorable 
geometries.  The  slope  of  (5.^6)  and  the  Joint  solution  of 
and  therefore  are  inci'eas Ingly  sensitive  to 

Instrument  and  computation  noise  as  the  maneuver  becomes 
small.  In  geometric  terms,  the  two  lines  in  the  (r,  r) 
space  cf  Figur’e  5.2  becoinc  par-allel  as  the  maneuver  becom.cs 
small . 

In  consequence,  it  is  evident  ti.at  the  s  i  adimet  ri  c/ li  near- 
kinematic  miethcd  outlir.ed  above  could  be  impreved  by  use  of  the 
nonlinear  kinematic  .method  described  in  Section  b. 
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b.  TECHNOLOGY  ASSESSMENT 

It  appears  that  linear  methods  of  kinematic  ranging 
can  work  effectively  only  for  favorable  combat  geometries 
and  if  the  sensor  noise  levels  are  quite  low.  This  Implies 
high  accuracy  sensors,  and  frequent  inef f ect Ivity . 

The  nonlinear  filter  described  in  Section  ^  estimates 
the  target  velocity-vector  components  and  therefore  can  be 
used  to  deduce  the  target  velocity  and  aspect  angle.  The 
filter  also  estimates  the  target's  acceleration  com.ponents, 
the  turn-rate, and  derivatjve.  It  therefore  estimates 

the  entire  specification  of  the  target  dynamics,  and  car,  be 
used  to  estimate  the  target  future  trajectory,  and  time  to  go. 
This  may  be  useful  not  only  to  the  aii'-air  missile  terminal- 
control  guidance  problem  but  also  to  a  variety  of  terminal 
fire  control  systems  in  aircraft  or  helicopters  instead  of 
missiles.  This  nonlinear  filter  can  also  be  used  to  eliminate 
impossible  target  accelerations  fromi  data  provided  by  other 
tracking  systems  such  as  radar  or  laser  rangers. 


Msny  of  th0  S06k.0rs  of'  TV  cr*  IR  typos  provide 

somie  Indication  of  the  target  sice  cr  shape.  Tai’get  size- 
data  can  be  used  as  a  partial  basis  for  a  ranging  system,  as 
outlined  in  Section  5.  In  addition.  Information  on  the 
target  shape  or  its  aspect-angle  can  be  used  to  provide 
clues  enabling  nonlinear  filtering.  For  example,  if  the 
target  is  viewed  from  a  head-on  or  tail-chase  aspect  it  is  clearly 

8h 
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difficult  fcr  the  target  to  accelerate  parallel  to  the 
Line  of  Sight;  it  can  accelerate  across  the  line  of  sight 
with  relative  ease. 

Linear  stadimetric  ranging  cannot  by  itself  estimate 
range  and  range-rate,  but  only  their  ratio.  However,  com¬ 
bing  linear  stadimetry  and  kinematic  schemeo  has  the  effect 
of  reducing  the  number  of  unknovm  variables  to  be  determined 
simultaneously.  The  parameter  identification  problem  may  in 
fact  be  reduced  to  a  problem  in  2-space,  without  using  the 

poor  assumptions  that  r,  a,p  ,  and  a„  are  quasi-constant. 

■^r  ^0 

In  this  situation,  tv;o  or  (preferably)  more  independent 
measurement  instants  enable  an  unique  solution.  The  problem 
of  updating  the  hyperplane  of  Figs.  5.1  or  5.2  to  another 
Instant  is  then  reduced.  It  is  not  eliminated,  therefore 
the  solution  is  not  rigorously  observable;  however,  the 
combination  may  be  useful. 

In  a  review  of  the  present  technology  and  an  estimate  of 
the  trends  and  potentials, the  Interaction  of  cost  with  systems 
technology  is  perhaps  most  important  in  the  air-air  missile 
problem.  These  points  .'tund  out; 

a)  Seekers  of  all  kinds  are  very  expensive;  a 
small  radar  seeker  is  very  costly  In  weight, 
space  and  dollars.  A  TV  seeker  may  be  less 
expensive  in  each  sense,  but  is  not  cheap.  The 
various  sophisticated  IR  seekers  are  cheaper  but 
are  still  quite  costly  compared  to  the  simplest 
IR  seekers. 

(b)  Sensors  such  as  rate  gyros  are  an  order 

of  magnitude  less  costly  than  simple  seekers. 
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c)  The  introduction  of  medium  or  large  scale 
integrated  digital  and  analog  chips  has  at 
once  Improved  the  reliability  of  the  compu¬ 
tation  process  and  reduced  the  cost  by 
orders  of  magnitude. 

Taken  together,  these  elements  suggest  that  cost/ef feet ivity 
considerations  lead  to  Improving  missile  system  performance 
by  combining  a  wide  range  of  the  inexpensive  sensors  v/itn  a 
relatively  sophisticated  computation  capability,  Thiis  may 
enable  use  of  a  relatively  unscphistlcated  seeker  of  rela¬ 
tively  low  cost. 

The  nonlinear  range/range  rate  estimator  shown  in 
Fig.  ^4.1  requires  a  considerable  computation  capability,  and 
3i^n3.ic  ^'r'crr 

also  requires  liie  missile  airspeed  components  which  may  be 
obtained  from  relatively  inexpensive  sensors  or  may  be 
estimated  from  accelerometer  and  rate-gyro  data.  But  the 
least  sophisticated  seeker  is  cufflcJent.  Our  approach  v;as 
to  maximize  the  information  which  could  be  extracted  from 
an  inexpensive  seeker  by  combining  several  inexpensive  sen¬ 
sors  with  a  relatively  sophisticated  yet  inexpensive 
computer  which  will  soon  be  available. 


v.n 
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APPENDICES 

APPENDIX  A 

The  Non linear  Estimator;  Derivation  and  Instrumentation 

It  Is  possible,  due  to  the  particular  foi'n.  of  the 
equations  and  constraints  of  the  Mathematical  Model  of  the 
system,  to  develop  a  nonlinear  I’ange  and  range-rate  identi¬ 
fier.  This  type  of  identifier  is  discussed  in  this  section. 

A . 1  Analysis 

The  key  elements  cf  the  Mathematical  Model  are  quoted 
from.  Section  3: 

Velocity  com.ponent  equations  parallel  to  the  Line-cf- 
Slght : 

f  -  cos  (y^-o)  +  Vj,^,  cos  (y  j,,-o )  =  0, 

r  -  +  Vj,  =  0. 

r  r 

(A.  1 

Velocity  component  equations  perpendicular  to  t.^.e 

LOS 

r  a  -  V^  sin  (Y,j,-a)  +  Vj,^^  sin  (Ym-o  )  =  Oj 
or  r  o  -  V„  +  V..  =0. 

1  i'l 
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We  also  demonstrated  that  the  acceleration  component 
equations  parallel  and  perpendicular  to  the  LOS  are, 


respectively : 


r  -  ro  -  a^  + 

r  r 


ro  +  2ro  -  +a„ 

a  a 


The  target  and  missile  velocity  and  acceleration 
components  parallel  and  perpendicular  to  the  LOS  are: 


Velocities : 


Parallel:  cos  i, 


Perpendicular:  sin 


(A.  3) 


(A. 


V  “V  sin  £ 
0 


Accelerations : 


(A. 5) 


Parallel:  a^p  =  cos  -  V^yt 


Perpendicular:  a^  =  sin  Cp  +  Vp  y-j.  cos 

o 


(A. 6) 


where 
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Angles : 


° 


Tm  -  ° 


(A 


The  principle  of  this  approach  in  sirtiply  tc  solve 
A , 3  for  r  as 

.  2 

r  =  ro  +  a,.  -  a^.  (A 

t"!  i 

r  r  • 

Integration  of  r  v.’ill  yield  i’  and  r  if  the  teriris  or.  th- 
rlrht  of  (A. 8)  are  available. 

The  term  ro  can  be  corupuied  easily  while  a..  ,  tie 

‘'r 

miss  lie  acceleration  parallel  to  tbio  LOG,  can  be  roasured 

directly.  The  problem,  thcrefcre  ir.  ncv;  to  o.stirriate  a.~ 

^  *  r 

in  some  v.-ay.  It  is  assurtied  that  initial  values  of  r  and 
are  provided. 

It  was  shown  that  the  target  acceleration  in  the  I'ange 

direction,  a.^.  ,  car.  ce  detcnnined  frcr.  available  data.  The 
^  r 

procedure  is  first  to  calculate  the  target  airspeed,  ty  u 

cf  (A.l)  and  (A. 2)  as 


( r  +  V, .  cos  4 

.'i 


;  +  tro 


V  J.  I . 


■■T 


cos  'C.n-  +  V,f  sin 


l  r. 


It  v.as  stated  ir.  Section  (m.1)  that  this  expression 


misy  be  di f f ei'ent iat ed  tc  yield  +  v_  a^.  =  , 

r  '  r  ‘  c  '  o  "  * 

This  is  den:on3trated  belov;.  Thei’efore,  let  us  differentiate 


(A.  9)  v;ith  respect  tc  tlm.e,  yielding 

2{(r+Vj,  cos  )  Lr-i- (\y  cos  sin  sir,  C,-] 

+  (ro+V._,  sin  C.^.;)[ro  +  fo  +(Vj/  sin  Cj.i+Vj.'f ,,  cos  t^,.)  - 


M  •> 
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Substitute  Tror.  (A.  5)  s  Imf  11  fy  it'g  notation  for  the 

acceleraticn  components  of  the  niissiie,  and  , 

''r  "c 

parallel  and  perpendicular  to  the  LOS,  respectively;  these 
accelerations  are  gathcrcc  in  parenthesis  v.-ithir.  the  saine 
brackets  in  the  expression  above  [(  )j. 


;  ^  o  j  i.  td  -1.  Ca  v' 


(i-  +  V,.  CCS  Cj-)(r  +  a,,  +  V,.o  sin 

r 

+  (ro  +  V,.  .  ir.  +  re  +  a,.  -  ’ 

1'!  !•  .1 

C 


Low  solve  (A.  2)  fer  V,.  sir.  i,,.,  niultif  ly  by  c  and  rei  Is-c*. 


»  i-  4-  ^  V.,.,  f  y  , 

U  i  i  ‘s-.  O  *w  A  1 1 1  \  •  r  • 


r  ^  •*'  1  ^  ^  T  '  V  .  •  V-,  , 


term  thus  fcrnied.  Lin.ilariy,  solV',  (A.l.  Tzi’  cos  , 
multiply  by  c  and  olitiiii.atc  the  teren  (V..  u  co?  C.,'  fro:;,  the- 
second  liiie  cf  (A.  10).  Those  subs  t  itucicnc  ;.’ield 

i.r  +  \1.  CCS  A.,,  ;(r  -  ro  +  a^.,  +cV:^  ni:.  a-;  + 


(re  +  V._.  Gin  Cj-}(rd  4  Arc  + 


-  0  V  _  c  ■; 


V  V  . 

'I  r  *  m  • 

;  X 


Keeping  the  rigne  side  ob  this  expi-essis'i! ,  gatiien  si 
coefficients  cf  V,j^.  on  the  lefi  ;  tfiis  '"ecuces  to 

U'+V  .  cos  Ar-,Mr  -ro  +  a,.  ;  +  r  ro  +  sin  t,..) 

J':  !m  .  .  I'j 

(rd  +2rc  +  a,.  ) {{r+V,.  ccs  £.-.)(V„  sin  a.t)  - 

i'l  !■:  J. 

G 


-(rc  +  V 


'j.-  sin  Cj.,)(V.,,  CCS 
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sin 

to  express  the  ta'T'get  velocity  componer.  3  parallel  and 
perpendicular  to  the  oOd . 

The  exact  solution  for  a^..  is  thus 

I 


r 
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Usually  Vrp  is  small  and  can  be  neglected,  and  all  other 
terms  on  the  right  of  (A. 15)  can  be  estimated. 

The  mechanization  of  this  algorithm  is  discussed  in 
Section  (A.  2),  follcv;ing. 

A.  2  Mechanization  of  the  Slg.ials  for  the 
_ Ideal  Algorithm _ 

The  range  acceleration  is  given  by  (A. 15);  it  may  be 
.  .  2 

integrates  to  yield  r  and  r.  The  term  ro  can  then  be 

formed,  v.'hile  the  missile  acceleration  component,  a^.  ,  can 

‘r 

be  directly  sensed. 

The  mechanization  of  (A, 15)  requires  that  r:;lssile 

airspeed  Vj_,,  the  missile  acceleration  components  aj,  and 

a,.  ,  and  the  angle  Cp,  be  available  as  signals  or  be  esti- 

o 

mated.  It  will  nov,'  be  shown  that  this  can  be  accomplished. 

( a )  Acceleration  Components 

The  acceleromieter  components  ,  a„  ,  and  a,,  ,  are  a 

i*l 

r  o 

special  case;  two  easy  solutions  are  available: 

(1)  Mount  the  accelerometers  on  tl-.e  seeker; 
then  one  yields  a.,  and  the  other  yields  a,. 

n  i'l 

r  o 

directly . 

(2)  Mount  the  accelerometers  on  the  missile 
body;  then  their  .signals  are  ap^  ,  and  ap.  ,  fore  - 
aft  and  cross-body,  respectively;  and  resolution 
of  the  signals  from  the  ac ce leronieter  through 


93 


Belt  Aerospace  Company 


the  "look-ancle,"  A,  yields 


“r' 

i  cos  A 

-sin 

A 

^sln  A 

cos 

(A. 16 


As  the  "look-angle,"  A,  between  the  missile 
cei'iterlirie  and  the  LOS  can  easily  be  r.,ad‘^r 
physically  available  in  a  seeker  gimbal ,  tiie 
resclution  of  the  acceleration  cemp ener.t n  is 
easily  accomplished. 


( b )  Velocity  Lst Imation 

If  the  rnlssile  body  I'itch  i-ate ,  c,  or  equivalent, 
j_ g n  v ?  1  ^ b  1 6*  s  1  T) ^  1 ,  ci n d  r u r t i r  ,  if  t  h c r* c c c  — 
body  acceleration  of  the  missile,  a.,  ,  is  also 

i: . . 

J 

available  as  a  signal,  it  is  relatively  easy  tc 
estimate  missile  velocity,  V._^. 

The  pertiiient  equations  of  inotlo.':  oi  i  he  mjlssile 


are  , 
a  = 


1*8 


l!u  +  G6 


v.'here  o  angle  of  ai.  Lack 

a^,  cross-tody  acceioratlon 

‘a; 

6  control  displacement 

q  body  pitch  rate 

V,.  airspeed 
1*. 
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i'J ,  G  cross-body  accelerations  f:er 

radian,  duo  to  ari(;:]e  of  attack 
aroj  control  disp lacercents  , 

I'espect.i  vely 

Assur':ji:f  cuasi-cor.otatil.  conditions,  tlio-  transl'ei' 
function  defiiiint^  acceleratior.  is 


a.. 

I'! 


s+:;/v,. 


T  r !  e  a  c  c  0 1  c  r  a  t  i  c  r.  d  u  o  t  o  c  c  n  t  r  o  1  dir-  p ! 
usually  to  nodUctoo  so  thoi 


iL 


I  U/X,.  \ 

‘  ‘V:  1  nrT:7v::l 


To  4  1 


(  A  .  1  c  } 


h.  •  j.  y  ; 


vvner’e 


Ti:-  !;  V  rriuy  bo-  ostiii.Cited  t,y  tt.o  a  1  y or-:;  t 


-  k 


/ 

V  U.  - 


I  • 


Q  ) 


or  V.,  -  o  (  |a-.  J  -  V..  |qi  : 

v;hore  V  >  0  has  th',-  ap;  rej  rlate  valu‘.  . 

oorr;ev.hat  suporior  ulrorilb!:, 
ar-o  aoriio'VOd  by  L-ja  -oinr  in':  alroilti,!.. 
lonyl t udinal  ac cclor'ut ion  c^mforioni  a,.  , 

I'l 

/. 


d.icii  vH.l]  lao  1- 


ivon  t  y  t;,o-  o.ioo:]' 


I,  .  1  CO 
(,'.0  1  ^ 
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Algorithm  (A. 21)  is  obviously  simple  to  instrument, 
requiring  only  one  multiplication. 

(c)  Estimation  of  the  Angle 

The  angle  C[,|  is  defined  as  • 

But  as  Yj,.  =  where  Bj,  is  the  centerline 

attitude,  ther.  4,;  =  (J,,-a-G ,  and  as  the 

angle  from,  m.issile  centerline  to  the  seeker  LOS  is 

X  =  0  -  6p 

we  have  4|j  = -X  -  a 

v.'here  a  is  the  angle  of  attack  betv/oen  the  n.issile 
centerline  and  the  velocity  vector.  But  this  can 
easily  bo-  estimated  as 


a 


Ts  +  1 


where  T  was  defined  above. 


(A.  23 


The  m.issile  velocity  components  and 

‘  r  ‘ 

can  then  b-  computed  as  Vj._.  =  cos  4.,,  tind 


=  V^,  Sin 
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APPENDIX  B 

Error  Analysis  Coefficients 

The  numerical  data  presented  in  this  appendix  represent 
the  time  varying  coefficients  of  the  error  equation.  As  dis¬ 
cussed  in  Chapter  4,  the  error  equation  is  of  the  form 

7 

X  +  a, X  +  a^x  =  I  b, N .  , 

and  the  data  here  represent  the  coefficients  and  tabulated 
for  every  0.25  seconds.  There  are  six  tables,  representing  a 
selection  of  six  choices  of  initial  geometry  configuration.  They 
are  ordered,  aivd  coi’respund  respectively  to  Figures  m.5  through 
^^.8  inclusive;  however,  each  graph  is  labelled  to  show  the  initial 
geometry.  These  data  permit  more  general  stochastic  analyses  of 
the  effects  on  the  estimates  due  to  errors  in  the  eral  signals 

which  are  used  in  the  algorithm.  The  coefficients  and  are: 

a^^  =  2a  tan  -  Y,p  tan  ,  and 
a^  =  -o  +  Trj-.c  +  o  tan 

The  coefficients  b^  and  erroi’  sources  are: 

^1  =  -1  >'3  »  a-celllLOS  error 

r 

=  sec  Crp  Vrj.,  target  airspeed  change 

b^  -Yrj.  sec  f;,p  sin  (Vj_.-Y.p)  Vj_, ,  missile  airspeed  error 
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^  = 

-rY,p  -2r 

tan 

+  2ro 

0 ,  LOS  rate  error 

^  = 

-Vt 

^T 

cos  (Yf.-Y.p) 

"5 

4]^;],  error  of  estimating  a 

-r  tan  C,j, 

0  LOS  acceleration  error 

b„  = 
1 

-tan 

a,.  ,  accel.  -1-  LOS,  error 

I'l  ’  ’ 

0 

All 

angles  are 

in 

radian  measure 

unles 

s  specifically  stated 

otherwise . 

The  data  of  this  appendix  are  presented  In  a  conventional 
"exponential  format"  for  digital  output  data.  Thus  the  tabulated 
item  0,^22¥^  01  is  interpreted  as  (0.^22)10'*'^  =  ^.22,  while  -0.2^^E-01 
is  interpreted  as  (-0.2*4^10  ^  =-0.02*4^4, 


H 

I 
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